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ABSTRACT 

NGC 7129 is a bright reflection nebula located in the molecular cloud complex near I =105?4, 
b =+9?9, about 1.15 kpc distant. Embedded within the reflection nebula is a young cluster 
dominated by a compact grouping of four early-type stars: BD-f65°1638 (B3V), BD-1-65°1637 
(B3e), SVS 13 (B5e), and LkHa 234 (B8e). About 80 Ha emission sources brighter than V ~23 
are identihed in the region, many of which are presumably T Tauri star members of the cluster. 
We also present deep (E~23), optical (VRcIc) photometry of a field centered on the reflection 
nebula and spectral types for more than 130 sources determined from low dispersion, optical 
spectroscopy. The narrow pre-main sequence evident in the color-magnitude diagram suggests 
that star formation was rapid and coeval. A median age of about 1.8 Myr is inferred for the 
Ha and literature-identified X-ray emission sources having established spectral types, using pre- 
main sequence evolutionary models. Our interpretation of the structure of the molecular cloud 
and the distribution of young stellar objects is that BD-|-65°1638 is primarily responsible for 
evacuating the blister-like cavity within the molecular cloud. LkHa 234 and several embedded 
sources evident in near infrared adaptive optics imaging have formed recently within the ridge 
of compressed molecular gas. The compact cluster of low-mass stars formed concurrently with 
the early-type members, concentrated within a central radius of ^0.7 pc. Star formation is 
simultaneously occurring in a semi-circular arc some ^3 pc in radius that outlines remaining 
dense regions of molecular gas. High dispersion, optical spectra are presented for BD-|-65°1638, 
BD-f65°f637, SVS 13, LkHa 234, and V350 Cep. These spectra are discussed in the context of 
the circumstellar environments inferred for these stars. 

Subject headings: stars: formation, pre-main sequence - Galaxy: open clusters and associations: individ- 
ual(NGC7129) 


1. Introduction 

The bright reflection nebula NGC 7129 {I = 
105°, 6 =-|-9?9) is illuminated by the early-type 
members of a young stellar cluster emerging from 
a molecular cloud complex composed of the dark 
nebulae Lynds 1181 (L1181), L1183, and TDS 
395 (Taylor et al. 1987; Yonekura et al. 1997). 
Signatures of star formation activity are present 
throughout the region including two dozen Herbig- 
Haro (HH) objects, two large molecular outflows, 
numerous nebulous filaments, and embedded in¬ 
frared sources. The cluster is dominated by 
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a tight grouping of B-type stars that includes 
BD-f65°1638 (B3), BD-b65°1637 (B3e), SVS 13 
(B5e), and the Herbig Be star LkHa 234 (B8e). 
An elongated cavity one parsec in diameter has 
been carved out of L1181, presumably by the stel¬ 
lar winds and ultraviolet flux emanating from the 
central grouping of B-stars. Shown in Figure 1 
is a three-color (Vi?cdc-band) image of the re¬ 
flection nebula obtained by G. H. Herbig in 1999 
October using the University of Hawaii (UH) 2.2 
m telescope on Mauna Kea. The image clearly 
depicts the dense molecular ridge to the east and 
south, the evacuated region around the four lumi¬ 
nous central stars, and the emerging young cluster 
lying within the reflection nebula. The distance of 
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NGC 7129 was estimated to be ~1 kpc by Racine 
(1968), but more recent estimates include 1.26 
kpc (Shevchenko & Yakubov 1989) and 1.15 kpc 
(Straizys et al. 2014), which is adopted in this 
analysis. Kun et al. (2008) present a comprehen¬ 
sive literature review of the Cepheus star forming 
region including NGC 7129 and its embedded clus¬ 
ter. Below we summarize only pertinent studies 
that are directly relevant to the present survey of 
the young stellar population. 

NGC 7129 and its associated molecular clouds 
have been studied extensively at millimeter wave¬ 
lengths by Loren (1975, 1977); Bechis et al. (1978); 
Bertout (1987); and Miskolczi et al. (2001). These 
^^CO and surveys have identified multi¬ 

ple density enhancements in the region, a semi¬ 
circular ridge of CO emission to the north, east, 
and south of the reflection nebula, and the large 
evacuated cavity centered near BD-|-65°1637. In 
the infrared, Bechis et al. (1978) detected two 100 
/rm sources in the area: LkHa 234 and far infrared 
source 2 (FIRS 2), which has no optical coun¬ 
terpart, located some 2' south of BD-|-65°1638. 
At least two molecular outflows have been iden¬ 
tified in the region (Edwards & Snell 1983; Liseau 
& Sandell 1983, Bechis et al. 1978); one cen¬ 
tered on FIRS 2 and the other originating near 
LkHa234. Yonekura et al. (1997) estimate the to¬ 
tal mass of the molecular cloud group associated 
with NGC 7129 to be ~2900 Mq. 

A single-channel 2.2 /rm survey of the region 
was carried out by Strom et al. (1976), here¬ 
after SVS, in a search for an embedded popula¬ 
tion of newly formed stars. Although none were 
identified, the detection limit of the survey was 
marginally close to the expected brightness of an 
average T Tauri star. Cohen & Schwartz (1983) 
identified four new infrared sources in the region, 
one of which was V350 Cep, an M2-type T Tauri 
star that brightened dramatically between 1970 
and 1976. Ibryamov et al. (2014) present a photo¬ 
metric study for this source that includes a long¬ 
term light curve constructed from archival data. 

Gutermuth et al. (2004) used the Two Micron 
All Sky Survey (2MASS), the FLAMINGOS spec¬ 
trograph on the Multiple Mirror Telescope, and 
the Infrared Array Camera (IRAC) onboard the 
Spitzer Space Telescope to complete an extensive 
1.2 to 8.0 /im survey of the cluster, ultimately de¬ 
termining a disk fraction of 54±14% within the 


cluster core (defined as a 0.5 pc radius arc ex¬ 
tending from the peak of the local stellar surface 
density). From a star count analysis, Gutermuth 
et al. (2004) estimate that there are about 80 clus¬ 
ter members within the core, down to a limiting 
magnitude of Ks=16 mag, or approximately 0.055 
Mq for unattenuated sources. In the peripheral 
areas Gutermuth et al. (2004) estimate that the 
extended cluster population is roughly equal to 
that of the core. 

Muzerolle et al. (2004) present 24 fim Multi¬ 
band Imaging Photometer for Spitzer (MIPS) 
imaging of NGC 7129 and provide a census of the 
cluster population outside of the reflection neb¬ 
ula. Of the 36 sources detected at 24 pm and in 
at least three IRAC passbands, 13 were classified 
as Class O/I protostars with envelopes including 
FIRS 2 (Class 0), another 18 as Class II sources 
(classical T Tauri stars, CTTS) with relatively low 
infrared luminosity, and six were classified as Class 
III sources with emission consistent with that of 
pure stellar photospheres. The majority of these 
sources are located in an arc sweeping north-east 
to south of the reflection nebula, effectively trac¬ 
ing the densest regions of the remnant molecular 
cloud. Muzerolle et al. (2004) conclude that star 
formation is not concentrated in the cluster core, 
but rather dispersed over a substantial area, some 
~3 pc in scale. 

Stelzer & Scholz (2009) combined the mid- 
infrared Spitzer surveys of Gutermuth et al. (2004) 
and Muzerolle et al. (2004) with a shallow 22 ks 
Chandra Advanced CCD Imaging Spectrometer 
(ACIS) X-ray observation of the region. From the 
Spitzer infrared excess sources and the 59 Chandra 
X-ray detections, they define a sample of 26 Class 
II and 25 Class III cluster members, estimated to 
be complete to ~0.5 Mq. The disk fraction of 
the least biased sub-sample composed of lightly 
extincted sources is '^33ljg%. 

While NGC 7129 has been studied extensively 
at millimeter, micron, and X-ray wavelengths, 
it remains poorly characterized in the optical. 
Hartigan & Lada (1985) obtained VRcIc and 
narrowband Ha CCD imaging of two regions in 
NGC 7129 which include most sources from the 
infrared surveys of Cohen & Schwartz (1983) and 
SVS. VRcIc photometry is provided for about 30 
sources in these two fields. Magakian et al. (2004) 
identified 22 Ha emission stars in the central and 
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northeast cluster regions, complete to V ~20 us¬ 
ing slitless spectroscopy. VRI— band photometry 
for about 100 sources was also obtained, but not 
provided in the published analysis. Magakian et 
al. (2004) concluded that most of the detected Ha 
emission sources were T Tauri stars. Vilnius 7- 
color photometry down to V =18.8 is provided by 
Straizys et al. (2014) for 159 stars in the cluster 
region. The photometric data was used to clas¬ 
sify about half of the detected sources, assign¬ 
ing spectral types, luminosity classes and deriv¬ 
ing extinction estimates. The extinction across a 
23'x23' area varies between 0.6 and 2.8 mag, but 
in the densest parts of the molecular cloud. Ay 
approaches ~13 mag. Straizys et al. (2014) deter¬ 
mine the age of the cluster to be between 0.2 and 3 
Myr by plotting six early-type cluster members on 
the HR diagram and using the evolutionary mod¬ 
els of Palla (2005). 

From the literature described above, some 90 
pre-main sequence candidates have been identi¬ 
fied in the cluster region spanning an area ~10' 
in diameter or ~3.3 pc in physical scale. These 
sources were detected using a number of youth 
discriminants including Ha emission (Hartigan & 
Lada 1985; Magakian et al. 2004), X-ray emis¬ 
sion (Stelzer & Scholz 2009), and infrared excess 
(Gutermuth et al. 2004; Muzerolle et al. 2004; 
Stelzer & Scholz 2009). No comprehensive optical 
photometric and spectroscopic census of the young 
cluster population, however, has been published to 
date. 

In this paper we present deep (V ~ 23) 
band photometry for a field roughly cen¬ 
tered on the reflection nebula. Results from a 
slitless grism, Ha spectroscopic survey of a ^^125 
square arcminute area of the cluster are discussed. 
Spectral types are provided for ~130 stars in the 
region, determined from low-dispersion, multi¬ 
object spectroscopy. We also present high disper¬ 
sion (R~45,000), optical spectra for BD-|-65°1638, 
BDH-65°1637, LkHa 234, SVS 13, and V350 Cep. 
Finally near-infrared, adaptive optics imaging 
of LkHa 234 is examined and compared with 
archived Spitzer IRAC data. 

This paper is organized as follows: in Section 
2 we discuss the observations made by the au¬ 
thors and by G. H. Herbig in support of this in¬ 
vestigation. In Section 3 we present the high- 
dispersion spectra of the early-type stars and of 


V350 Cep; in Section 4 the stellar population of 
NGC 7129 is examined to include the early-type 
stars, the Ha emission sources, and the X-ray de¬ 
tected sources. We present in tabular form pho¬ 
tometry and spectral types determined at optical 
wavelengths for two subsets of stars: 1) those ex¬ 
hibiting Ha and/or X-ray emission, and 2) others 
in the field, some of which exhibit infrared excess. 
In Section 5 we discuss the color-magnitude di¬ 
agram of the pre-main sequence population, the 
average extinction suffered by those stars with as¬ 
signed spectral types, the ages and masses of the 
cluster members, and the infrared color-color dia¬ 
grams. In Section 6 we discuss the progression of 
star formation in the parent molecular clouds as 
well as the star forming history in the region. We 
examine the source of the molecular outflow near 
LkHa 234 and the origin of the photo-dissociation 
region enveloping the early-type stars. 

2. Observations 

2.1. Optical Photometry 

Three epochs of optical imaging data were ob¬ 
tained of the NGC 7129 region between 1993 and 
1999. Extinction corrections and transformation 
to the Landolt (1992) system were achieved for 
all three sets by observing several standard fields 
periodically throughout the night at various air 
masses. The photometric calibrations included 
zero point, airmass, and color terms except when 
indicated below. 

The first and most extensive survey consisted 
of V Rc^C—band CCD imaging obtained using the 
T2KA CCD at the Kitt Peak National Observa¬ 
tory (KPNO) 0.9m telescope in 1993 June by LAH 
and S. Strom. The T2KA detector is a 2048x2048 
CCD with 24 /xm pixels yielding a platescale of 
0769 pixel“^. Twilight flats were obtained at the 
beginning and end of the nights in all filters to 
remove pixel to pixel variations in response. Ex¬ 
posure times were 3, 30 and 300 s in all passbands. 
The optical images span a region 23'x23' in area 
with a photometric completeness limit of V ~18.5. 
The faintest sources evident on the CCD images 
are about three magnitudes fainter. The region 
imaged by the KPNO T2KA CCD encompasses 
the 20'x20' red Digitized Sky Survey (DSS) im¬ 
age of the region shown in Figure 2. The T2KA 
CCD images were reduced using standard tasks 
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and procedures available in the Image Reduction 
and Analysis Facility (IRAF). Point spread func¬ 
tion (PSF) fitting photometry was performed us¬ 
ing the DAOPHOT package of IRAF. The V — Ic, 
V color-magnitude diagram for all ~2500 sources 
with photometry available from the T2KA CCD 
imaging is shown in Figure 3 (left panel). We pro¬ 
vide a supplemental on-line catalog with the J2000 
coordinates, optical (V, V — Rc, and Rc — Ic) 
photometry, and photometric errors for all ~2500 
sources. 

The second epoch of band photometry 

was obtained using the Low Resolution Imaging 
Spectrograph (LRIS) on Keck II in 1999 June by 
LAH. The 8' x6' LRIS images were centered on 
LkHa 234 (Figure 2) and reveal dozens of faint 
cluster members within the compact grouping of 
early-type stars. Exposure times ranged from 3 to 
600 s to ensure that the majority of sources were 
unsaturated. Tasks within the PHOT package of 
IRAF were used to carry out aperture photometry. 
Photometric calibration was achieved using zero 
point and airmass terms only, i.e. no color terms 
were applied. Shown in Figure 3 (right panel) is 
the Rc — Ic, Rc color-magnitude diagram for all 
sources in the LRIS field of view with no allowance 
being made for interstellar reddening. 

Finally BVRcIc photometry of the NGC 7129 
region was obtained in 1999 October at the f/10 
focus of the UH 2.2 m telescope on Mauna Kea 
by G. H. Herbig. The detector was a Tektronix 
2048x2048 CCD with 24 /xm pixels. The field was 
imaged in all filters during photometric conditions 
and in average seeing conditions for Mauna Kea, 
0"75. The plate scale was 0'.'22 pixel“^, and the 
imaged field was approximately 7.'5x7.'5 in area. 
Three exposures were obtained in each filter with 
integration times of 5, 60, and 300 s. Image reduc¬ 
tion and analysis were accomplished by SED using 
IRAF. Aperture and PSF fitting photometry were 
carried out with the PHOT, PSF, and ALLSTAR 
tasks in the DAOPHOT package of IRAF. Bright, 
isolated, single stars were selected across the field 
to construct a characteristic PSF for the images 
in each color. 

Photometry from the three optical imaging 
surveys are compared in Figure 4, 1993 June 
(KPNO) vs 1999 October (UH 2.2 m) and Figure 
5, 1993 June (KPNO) vs 1999 June (Keck LRIS). 
There is relative agreement, however a slight offset 


(~-|-0.16 mag) exists between the V — Ic colors 
of the 1999 October UH 2.2 m observations and 
the 1993 June KPNO observations. This offset 
originates from extended red transmission from 
the UH 2.2 m Ic filter, which is documented by 
Courtois et al. (2011). Given the presence of 
this red leak in the UH 2.2 m Ic filter, the optical 
photometry presented here is taken primarily from 
the 1993 KPNO observations, unless photometry 
for a given source were not available. In these 
circumstances, the photometry from the 1999 UH 
2.2 m observations were substituted with a —0.16 
mag correction being applied to the V — Ic color. 
Differences in the Rc-Ic colors of the 1993 KPNO 
and the 1999 Keck LRIS data evident in Figure 
5 likely arise from the missing color term in the 
photometric calibration of the latter data set. 

2.2. Ha Slitless Grism Spectroscopy 

The Ha emission survey of NGC 7129 was car¬ 
ried out on 10 October 2003 by SED using the 
Wide-Field Grism Spectrograph (WEGS) installed 
at the f/10 Cassegrain focus of the University of 
Hawaii 2.2 m telescope. A 420 line mm“^ grism 
blazed at 6400A yielded a dispersion of 3.85A 
pixeU^. The narrowband Ha filter isolated a 
region of the first-order spectra between ~6300 
and 6750A. The WEGS imaged upon the central 
1024x1024 pixels of a Tektronix 2048x2048 pixel 
GGD, yielding a nominal field of view ~5f5x5f5. 
The WEGS survey covered a region of about 150 
square arcminutes, roughly centered upon LkHa 
234. Integration times for the slitless grism spectra 
were 120 and 1200 s for each field to ensure com¬ 
plete, unsaturated coverage for nearly all sources. 
The minimum measurable equivalent width of the 
WFGS spectra is about 2 A. 

2.3. Low Dispersion Optical Spectroscopy 

Low-dispersion spectra were obtained using 
HYDRA on WIYN in 1993 and the RC Spec¬ 
trograph on the Mayall 4 m telescope in 1994 by 
LAH. Spectra of earlier-type stars (<K5) were 
classified using the standards of Allen & Strom 
(1995). The relative strengths of a number of 
temperature sensitive features were adopted to in¬ 
clude the Na I D doublet, the Ca I absorption lines 
between 6100 and 6200A, the blend near Ha, and 
the Ca II near-infrared triplet, if in absorption. 
For late-type stars (K5 and later), the standards 
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of Kirkpatrick et al. (1993) and the temperature 
sensitive TiO indices of Slesnick et al. (2006) and 
references therein were adopted. These indices 
measure the strength of TiO absorption features 
near A7140 and A8465 relative to narrow contin¬ 
uum passbands at 7035 A and 8415 A respec¬ 
tively. The resulting indices were then compared 
with those of dwarf and pre-main sequence objects 
compiled by Slesnick et al. (2006). Shown in Fig¬ 
ure 6 are examples of the WIYN spectra, ordered 
by spectral type with the regions used by the TiO 
indices for classihcation purposes demarcated in 
gray. The spectra were independently classified by 
both authors with LAH following methods similar 
to those described in Hillenbrand (1997). Where 
type or specific subclass could not be determined, 
ranges were estimated and are provided in Tables 
1 and 3. 

2.4. High Resolution Echelle Spectrome¬ 
ter (HIRES) Observations 

High-dispersion optical spectra were obtained 
by G. H. Herbig for BD-h65°1638, BD-H65°1637, 
SVS 13, and V350 Cep using the High Resolu¬ 
tion Echelle Spectrometer (HIRES; Vogt et al. 
1994) on the Keck I telescope during the nights 
of 06 July 2003, 24 July 2004, 06 June 2010, and 
22 July 2007, respectively. The HIRES spectra 
of LkHa 234 were obtained by LAH on 06 Dec 
1999 and by J. Kuhn on 13 Jun 2004. HIRES 
was configured with the red cross-disperser and 
collimator in beam for all observations presented 
here. The Cl decker (0'.'87x7'.'0), which has a pro¬ 
jected slit width of 3 pixels, was used for the bulk 
of the spectra, providing a spectral resolution of 
~45,000 (~6.7 km s“^). Near-complete spectral 
coverage from ~4340 to 6800A was achieved for 
the 2003 and 2004 observations, which were made 
with the single Tektronix 2048x2048 pixel CCD 
camera with 24 y^m pixels. The 1999 spectrum of 
LkHa 234 used the D2 decker (I'.'148x2870) and 
covers a range from 6240 to 8680 A. The 2007 
observations of V350 Cep were made with the up¬ 
graded three-chip mosaic of MIT-LL CCDs having 
15 fj,m pixels. The resulting spectral coverage of 
the 2007 observations of V350 Cep range from ap¬ 
proximately 4300 to 8600 A. The CCDs were used 
in low-gain mode, resulting in readout noise lev¬ 
els of ~2.8, 3.1, and 3.1 e“^ for the red, green, 
and blue detectors, respectively. Internal quartz 


lamps were used for flat fielding and ThAr lamp 
spectra were obtained for wavelength calibration. 
Integration times were 120-1200s, yielding signal- 
to-noise levels of up to '^lOO for the early-type 
stars. The cross-dispersed spectra were reduced 
and extracted using the Mauna Kea Echelle Ex¬ 
traction (MAKEE) pipeline written by T. Bar- 
low. Standard routines available through IRAF 
and IDL were used for spectral analysis. The 
raw as well as pipeline extracted spectra for the 
HIRES observations presented here are available 
through the Keck Observatory Archive (KOA): 
https://koa.ipac.caltech.edu. 

2.5. Adaptive Optics Imaging of LkHa 234 

High angular resolution, near infrared imaging 
of LkHa 234 was obtained with NIRC2 using the 
natural guide star (NGS) adaptive optics system 
on the Keck H telescope on 2012 June 30 by SED. 
The images were obtained with the wide cam¬ 
era (0704 pixel“^ platescale) in beam yielding a 
40"x40" field of view. A three position dither 
was used to image the field, avoiding the quadrant 
of the Alladin HI detector that exhibits a higher 
level of fixed pattern noise. Integration times were 
10 s with 10 coadds in all three filters {J,H,K'), 
yielding an effective integration time of 100 s per 
frame. Basic image reduction and analysis were 
completed using standard routines written in IDL. 

3. Analysis of the High Dispersion Spec¬ 
troscopy in the Context of Circumstel- 
lar Environments 

3.1. BD-|-65°1637 

BD-|-65°1637 was recognized by Merrill & Bur- 
well (1950) as a moderate intensity Ha emission 
star. Herbig (1960) found similar emission on 
plates obtained at Lick Observatory and classi¬ 
fied the star as B5n or earlier. Narrow emission 
components were present near the centers of all 
Balmer lines from H/3 through HJ and possibly 
extending to He. Herbig (1960) also noted that 
interstellar Ca II was present and that the stellar 
spectrum was suggestive of an ordinary Be star. 
Other classifications for BD-|-65°1637 taken from 
the literature include: B2-B3 (Strom et al. 1972), 
B5 (Einkenzeller 1985), and B3 (Hillenbrand et 
al. 1992). Strom et al. (1972) found H/3 to be in 
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emission with no obvious P Cygni component and 
noted that BD+65°1637 is a very rapid rotator. 

Using the He I A4471/Mg II A4481 line ratio 
found in the HIRES spectrogram of BD+65°I637, 
we assess the temperature class to be ~B3, con¬ 
sistent with classifications found in the literature. 
The HIRES spectrum reveals strong Ha emis¬ 
sion (IU=—25 A) with wings extending out to at 
least ±550 km s“^. There is no indication of an 
underlying, early-type stellar photosphere in the 
Ha emission profile. The emission line is double- 
peaked as shown in Figure 7 with radial veloci¬ 
ties measured for the blue-shifted and red-shifted 
components of —60.5 and ±13.7 km s“^, respec¬ 
tively. The radial velocity of the core absorption 
feature is —29.8 km s“^. These velocities differ 
substantially from those found by Finkenzeller & 
Jankovics (1984), —5 km s“^ (core absorption) 
and —36 and ±45 km s“^ (blue and red emission 
peaks, respectively). H/3, shown in Figure 7, re¬ 
veals a similar double-peaked emission structure, 
but is substantially weaker (1U=—2.4 A). The ve¬ 
locities of the double-peaked H/3 emission profile 
are —78.6 and ±30.1 km s“^, with the red emis¬ 
sion component having a higher amplitude than 
the blue component. Unlike the Ha profile, broad 
absorption wings are present in H/3. Hy falls out¬ 
side the coverage of the HIRES spectrum, but the 
wing of its redward absorption edge is evident on 
the bluest order. If Hy were in emission, it would 
be expected to be weak. 

Other features of interest in the HIRES spec¬ 
trum of BD±65°1637 include P Cygni-like pro¬ 
files of He I A4922 and A5015 (possibly contam¬ 
inated by Fe II emission), and broad, double- 
peaked emission profiles of Fe II AA5169, 5316, 
6317 and 6456. These iron features consistently 
exhibit higher amplitude red emission peaks rela¬ 
tive to blue (similar to Ha and H/3). The Na I D 
lines (AA5889.95, 5895.92) exhibit no trace of emis¬ 
sion and the narrow interstellar absorption compo¬ 
nents have heliocentric radial velocities of 19.7 
km s“^. In the local standard of rest, these veloci¬ 
ties are consistent with being produced within the 
local spiral arm, no more than 1 kpc distant using 
the galactic rotation curves calculated by Miinch 
(1957). Diffuse interstellar bands (DIBs) are also 
evident in the spectrum near AA5780, 5797, 5825, 
and 5850. 

The heliocentric radial velocity of BD±65°1637 


is difficult to measure given the very broad ab¬ 
sorption profiles induced by rapid rotation and the 
presence of Fe II emission near select He I absorp¬ 
tion features. Alecian et al. (2013) determined 
a usim for BD±65°1637 of 278±27 km s“^ and 
a radial velocity of — 26±20 km s“^, consistent 
with the heliocentric radial velocity of the molec¬ 
ular cloud reported by Finkenzeller & Jankovics 
(1984), —23.5 km s“^. Using a handful of He I 
absorption lines, we estimate the radial velocity of 
BD±65°1637 to be -25.6 km s"!. 

Hillenbrand et al. (1992) classify BD±65°1637 
as a Group HI object, i.e. a Be star with substan¬ 
tially less infrared excess than sources associated 
with massive envelopes or circumstellar disks. The 
infrared excess exhibited by BD±65°1637 likely 
results from free-free emission from hot circum¬ 
stellar gas, comparable to classical Be stars. Its 
association with the star forming region, however, 
implies that the star is young. In their mid and 
far-infrared Infrared Space Observatory (ISO) sur¬ 
vey of Herbig AeBe stars, Abraham et al. (2000) 
found that BD±65°1637 was weakly detected at 
60 /im, contrary to expectations if indeed disk- 
free. If conhrmed, Abraham et al. (2000) suggest 
that this 60 ^m emission could arise from cold 
dust associated with the star. Lorenzetti et al. 
(2003) mapped NGC 7129 using the Long Wave 
Spectrometer onboard ISO in [O I] 63 and 145 /rm 
and in [G H] 158 /rm. The line emission is sugges¬ 
tive of two photodissociation regions (PDR), one 
being illuminated by BD±65°1637 and the other 
by LkHa 234. In their Spitzer MIPS 24 fim sur¬ 
vey of the region, Muzerolle et al. (2004) report 
that while the photospheres of BD±65°1637 and 
BD±65°1638 were not detected directly, extended 
halos of warm dust are apparent around their po¬ 
sitions. The MIPS observations appear to conhrm 
that a halo of dust enshrouds these early-type stars 
on scales larger than an envelope / disk environ¬ 
ment. 

3.2. BD±65°1638 

BD±65°1638, was classified by Hubble (1922) 
as B3, and is among the earliest spectral type 
cluster members. Racine (1968) also assigned a 
B3 type and derived an absolute magnitude of 
My = —2.0, somewhat more luminous than a nor¬ 
mal B3 zero age main sequence star. The slightly 
greater distance adopted here would imply that 
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the star is even more luminous. 

Matthews et al. (2003) note that BD+65°1638 
is accompanied by a small H II region and that 
its temperature and luminosity place it near the 
birthline of a ~6 Mq star using the evolutionary 
models of Palla & Stabler (1993). Matthews et al. 
(2003) assign a spectral type of B2.5 based upon 
1.4GHz continuum flux and the star’s estimated 
excitation parameter. 

The HIRES spectrogram of BD+65°1638 shows 
a He I A4471/Mg II A4481 line ratio that is con¬ 
sistent with a B3 spectral type. The spectrum 
also reveals possible evidence for emission reversal 
within the core of the Ha absorption profile shown 
in Figure 8. The centroid of this weak feature has 
a velocity of about —25.8 km s“^, which is con¬ 
sistent with the radial velocity of the molecular 
cloud. H/3 appears to be a normal absorption pro¬ 
file (Figure 8) as might be expected for an early- 
type stellar photosphere. No emission features at¬ 
tributed to Fe II are present in its spectrum as 
with BD-|-65°1637 and LkHa234 (discussed be¬ 
low). 

The measured radial velocities for the interstel¬ 
lar Na I D lines are —21.9 and —21.3 km s“^; 
comparable to those measured for BD-|-65°1637. 
The heliocentric radial velocity of the star mea¬ 
sured from Mg II A4481, He I AA4387, 4437, 4471, 
4713, 4921, 5015, 5575, 5875 and 6678 is estimated 
to be ~-|-8.1 km s“^. Given the large disparity in 
the radial velocity of BD-|-65°1638 relative to the 
molecular cloud, it is likely that BD-|-65°1638 is 
a spectroscopic binary. Follow-up high dispersion 
spectroscopy is needed for confirmation. Other¬ 
wise the optical spectrum of BD-|-65°1638 is quite 
unremarkable, particularly if the star were lying 
on or near the stellar birthline. 

3.3. LkHa 234 

The luminous Herbig Be star LkHa 234 was 
classified by Herbig (1960) as having a late-A spec¬ 
tral type with strong Ha emission and moderate 
intensity H/3 emission. Strom et al. (1972) clas¬ 
sify the star as B5-B7 and noted that H/3 and Hy 
were in emission with sharp P Gygni-like com¬ 
ponents. Possible weak emission was also noted 
for He I AA4922, 5015, and the star was identi¬ 
fied as a moderately slow rotator. More recently 
Hernandez et al. (2004) assign a spectral type 


of B7 for LkHa 234. The HIRES spectrum of the 
source obtained in 2004 reveals a He I A4471 /Mg II 
A4481 line ratio that is consistent with a B8 spec¬ 
tral type. 

LkHa 234 is embedded within the ridge of 
molecular gas that defines the northeastern edge 
of the evacuated cavity. From its placement above 
the zero age main sequence in the color-magnitude 
diagram, it is generally assumed that LkHa 234 
is among the youngest of the massive stars in 
NGG 7129. Hillenbrand et al. (1992) classify it as 
a Group I object with substantial infrared excess 
beginning near 1 /rm and extending out beyond 
100 /rm. 

The HIRES spectra of LkHa 234 obtained in 
1999 and 2004 are suggestive of active accretion 
with numerous metallic and forbidden transitions 
in emission as well as Ha and H/3 exhibiting com¬ 
plex emission structure. Shown in Figure 9 are 
the Ha profiles for LkHa 234 from 1999 and 2004. 
While of comparable strength (LF=—56.6A), the 
profiles are significantly different in appearance 
with the 1999 spectrum exhibiting deep absorp¬ 
tion in a P Cygni-like profile. A second narrow 
absorption core is present near line center, just 
blueward of the sharp-edged emission peak. This 
peak and its stepped redward slope are remarkably 
similar in the two observations. The Ha profiles 
possess broad wings that extend several hundred 
km s“^ from the central wavelength. The blue- 
ward side of the 2004 emission peak shows a dis¬ 
tinct, sharp drop-off that never declines below the 
continuum level (Figure 9). The Ha profile is sug¬ 
gestive of a wind moving outward from the central 
star. This is supported by strong, blue-shifted [O 
I] AA6300, 6363 emission with heliocentric radial 
velocities of abont —23 km s“^, as well as blue- 
shifted [S H] AA6717, 6731 emission with compa¬ 
rable velocities. The radial velocity of LkHa 234 
could not be determined from the HIRES spec¬ 
tra as individual He I and Mg II absorption lines 
yield inconsistent results, possibly resulting from 
significant stellar activity. 

Other prominent features in the spectrum 
of LkHa 234 include complex profiles of Fe II 
AA4922,5018, which appear to possess broad emis¬ 
sion bases with central absorption and core emis¬ 
sion reversal rising well above the continuum level. 
The emission core has a radial velocity approxi¬ 
mately equal to that of the other forbidden emis- 
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sion lines, implying that it is associated with an 
outward moving wind. The Na I D lines exhibit 
narrow and broad absorption components, the 
former with radial velocities of about —20 km s“^ 
(interstellar) and the latter, —83 km s“^. Shown 
in Figure 10 are the profiles of H/3, Fe II A5018, 
and the Na I D lines in the 2004 HIRES spectrum 
of LkHa 234. The radial velocities associated with 
various features in the spectrum are annotated for 
reference. 

3.4. SVS 13 

SVS 13 was previously unclassified in the liter¬ 
ature, but the high dispersion spectroscopy pre¬ 
sented here is suggestive of a mid- to-late B spec¬ 
tral type with very broad He I AA492I, 5875, 6678 
absorption lines, implying rapid rotation. The 
profiles of H/3 and Ha are shown in Figure 11. 
While H/3 appears to be in pure absorption, Ha 
reveals a broad absorption profile with a narrow 
emission core having a heliocentric radial velocity 
of —23 km s“^. Fe I AA4988, 4994 are in weak 
emission as is a broad [O I] A6300 feature. The 
Na I D lines exhibit no trace of emission and the 
narrow interstellar absorption components have 
heliocentric radial velocities of 19.5 km s“^, 
consistent with those found in the other early-type 
stars in the cluster. Diffuse interstellar bands are 
also present near AA5780, 5797, and 5850. 

Stelzer & Scholz (2009) identify SVS 13 as an 
X-ray source (S3-X51) and from its 2MASS and 
IRAC [3.6], [4.5] photometry, as a Class HI in¬ 
frared source, i.e. consistent with a stellar photo¬ 
sphere. Muzerolle et al. (2004), however, find the 
PSF core of SVS 13 in their MIPS 24/im imag¬ 
ing to be completely saturated. The AT-band 
photometry presented here was obtained with the 
SQHD camera on the KPNO 50-inch telescope in 
1993 (2MASS only provides a Ks upper limit for 
this source). The star’s placement in the near- 
infrared color-color diagrams to be discussed in 
§5.4 is consistent with circumstellar dust emission. 
If the observed infrared excess arises from a cir¬ 
cumstellar disk, the Ha emission evident in the 
HIRES spectrum would be consistent with that of 
a Herbig Be star. 


3.5. V350 Cep 

The late-type (M2) pre-main sequence star 
V350 Cep lies near the edge of the L1181 molecu¬ 
lar cloud, ~5l0 north of the cluster core. Around 
it lie several Ha emission sources as well as the 
bright infrared source SVS 10. Herbig (2008) pro¬ 
vides an abridged history of V350 Cep observa¬ 
tions, which include that the source was unde¬ 
tected on the 1954 Palomar Sky Survey plates, 
but rose to its current brightness (V~16) some¬ 
time in the late 1960s or early 1970s. The source 
is clearly visible on plate 1 of SVS, a red photo¬ 
graph of NGC 7129 obtained using the Mayall 4 
m telescope at Kitt Peak. Pogosyants (1991) sug¬ 
gests that V350 Cep first exceeded the brightness 
limit {V ^18.5) of archived photographic plates 
sometime during 1971. Following its apparent 
minimum, V350 Cep has now remained “bright” 
for over 40 years, and it remains unclear if this 
activity is repetitive as might be expected of an 
EXor. 

Spectra of V350 Cep obtained by Magakyan & 
Amirkhanyan (1979) suggest a T Tauri-like source, 
but were clearly variable as described by Maga¬ 
kyan (1983). Cohen & Fuller (1985) found strong 
Balmer line emission as well as [O I] and Fe II emis¬ 
sion with an underlying M2-type photosphere, in¬ 
ferred from shallow absorption features in the red 
thought to arise from TiO bandheads. Hartigan & 
Lada (1985) conclude that V350 Cep illuminates 
the small reflection nebula GGD 33 and may be an 
outflow source within the region. Goodrich (1986) 
obtained a low-resolution, spectrophotometric ob¬ 
servation of V350 Gep and noted strong Balmer 
line emission as well as a rich Fe II emission spec¬ 
trum. 

The HIRES spectrum of V350 Cep presented 
here is described briefly by Herbig (2008), who 
otherwise focuses his discussion on an earlier 2005 
spectrum, as exhibiting strong Balmer line emis¬ 
sion as well as metallic line emission. The promi¬ 
nent metal lines are evident throughout the spec¬ 
trum including Fe I AA6I36, 6191, 6230 (Figure 
12); Fe II AA6317, 6338, 6456. Strong He I emis¬ 
sion is present including AA5875, 6678. The Na I 
D lines (Figure 12) exhibit broad emission with 
narrow interstellar absorption components within 
their cores having radial velocities of —18 km s“^, 
comparable to the velocities measured in the early- 



type stars. Broad Na I absorption features are also 
evident shifted some —204 km s“^ blueward of the 
rest wavelength. Herbig (2008) suggests that these 
features originate from the same outflowing ma¬ 
terial responsible for the P Cygni feature in Ha. 
Similar radial velocities were reported by Herbig 
(2008) for these broadened absorption features. 

Forbidden emission, e.g. [O I] A6300, [N HI] 
A6515, [N H] A6545, [O H] A6589, is present 
throughout the HIRES spectrum. Lithium A6708 
is in absorption with an equivalent width of 
~0.10A, slightly reduced from that measured by 
Herbig (2008), ~0.I4A in the earlier HIRES spec¬ 
trum obtained in 2005. The heliocentric radial 
velocity of the star determined using a dozen 
metallic features in the 5700-6800A wavelength 
region is —20.7±0.9 km s“^, consistent with that 
measured by Herbig (1998) using the 2005 HIRES 
spectrum, —23.0±0.3 km s“^. 

Ha exhibits a strong P Cygni-like profile, W = 
—27A with a weak absorption feature lying within 
the core of the emission peak that is not evident 
in the 2005 spectrum obtained by Herbig (2008). 
The deep P Cygni absorption structure at Ha has 
edges of —102 and —252 km s^, defined at the 
continuum level. Figure 12 shows the Ha emis¬ 
sion profile with the velocities annotated for vari¬ 
ous features. This figure can be directly compared 
with that shown in Figure 9 of Herbig (2008) from 
the 2005 HIRES spectrum. The overall shape of 
the line profile has changed somewhat while the 
equivalent width decreased only slightly. 

Herbig (2008) concludes that V350 Cep is not 
an EXor candidate, i.e. a T Tauri-like star that 
undergoes periodic flare-ups presumably as the re¬ 
sult of an inflow of material from the surrounding 
accretion disk. Muzerolle et al. (2004) classified 
V350 Cep as a Class H source from its placement 
in the Spitzer [3.6] —[5.8], [8] —[24] color-color di¬ 
agram. While eliminated as an EXor candidate, 
V350 Cep is clearly undergoing significant accre¬ 
tion activity and is associated with a substantial 
circumstellar disk as inferred from its strong Ha 
emission and infrared excess. The nature of it dra¬ 
matic rise from obscurity is worthy of additional 
study. 


4. NGC7129: The Cluster Population 

4.1. The Early-Type Stars 

The four central, early-type stars BD-|-65°1638, 
BD-f65°1637, LkHa 234, and SVS 13 are the most 
prominent members of the emerging young cluster 
in NGC 7129. In addition to these sources, how¬ 
ever, are approximately nine other bright stars in 
the region, dispersed over several arc-minutes on 
sky. 

On the southwestern edge of the evacuated cav¬ 
ity lie BD-t65°1635 (SVS 5) and BD+65°1636. 
The former is a G8-type star with measured proper 
motion and a probable field interloper, while the 
latter is a B8-type X-ray source presumably as¬ 
sociated with the young cluster. Six arcminutes 
northwest of BD-|-65°1638 lies BD-|-65°1631, aKO- 
type giant lying in the foreground. Along a sim¬ 
ilar axis, 3f5 distant from BD-|-65°1638, lies SVS 
4, a field G-type giant. Several arcminutes north 
of the cluster, three moderately bright nebulous 
stars form an arc terminating near V350 Cep: SVS 
2 (Al), SVS 10 (B8) and SVS 11. All of these 
sources are probable pre-main sequence stars that 
appear to be involved with the nebulosity. To the 
northeast, about six arcminutes from the cluster 
core, lies 2MASS J21435035-f6608477, classified 
by Straizys et al. (2014) as a B7 dwarf. This source 
appears to illuminates a small reflection nebula on 
the edge of the molecular cloud. Finally, four ar¬ 
cminutes southeast of BD-|-65°1638 is SVS 16, a 
G6-type star and X-ray source that exhibits an un¬ 
remarkable absorption line spectrum. This source 
appears to lie in the foreground of the molecular 
cloud, and its membership in the cluster remains 
unconfirmed. 

In addition to these, the intermediate mass 
(~2-8 Mq) Class 0 source FIRS 2, is still embed¬ 
ded within its natal envelope south of the cluster 
core and will likely emerge as a mid-B type star 
(Fuente et al. 2014). In summary, there are at 
least four other B-type stars that are associated 
with NGC 7129 and the molecular cloud complex. 

4.2. The Ha Emission Sources 

In his landmark study of Be and Ae stars as¬ 
sociated with nebulosity, Herbig (1960) reported 
that both BD-|-65°1637 and LkHa 234 exhibit Ha 
emission, but also noted the presence of several 
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faint stars having Ha emission near the limit of 
the slitless grating exposures within the nebulos¬ 
ity of NGC7129. From R/Ha magnitude ratios, 
Hartigan & Lada (1985) found that in addition 
to LkHa 234 and BD+65°1637, their sources 19S, 
14S, 26S, 13N (V350 Cep), 19N, and 27N were 
possible emitters. Several HH objects were also 
either confirmed or identified from the Ha imag¬ 
ing including GGD 32, 34 and 35, HH 103 and HH 
105. Hartigan & Lada (1985) suggested that LkHa 
234 is the driving source for at least some of these 
HH objects. Magakian et al. (2004) identified 22 
emission-line sources in NGC 7129, 16 of which 
were previously unknown. G. H. Herbig provided 
Magakian et al. (2004) with a list of Ha emission 
sources identified on an earlier WFGS image (not 
available for this analysis). Of 13 sources identified 
by Herbig as exhibiting Ha emission, five were not 
found in emission by Magakian et al. (2004), while 
their source MMN 11 was not identified by Her¬ 
big as an emitter. The WFGS images examined 
here reveal what appears to be a flat continuum 
for MMN 11 (i.e. no absorption or emission), con¬ 
sistent with Herbig’s findings. Of the 22 emission¬ 
line stars, Magakian et al. (2004) found that about 
half were concentrated in the central region of the 
reflection nebula. Ten of their sources were con¬ 
fidently identified as CTTS and another seven as 
weak-line T Tauri stars (WTTS). 

The Ha slitless grism survey presented here 
identihed over 50 Ha emission sources in a region 
some 150 square arcminutes in area, which is out¬ 
lined in Figure 2. The HYDRA low-dispersion 
spectroscopy revealed ~30 additional emitters, 
most of which lie outside the boundaries of the 
WFGS survey, have measured equivalent widths 
below the detection threshold of the WFGS (~2 
A), or had weak continua on the WFGS images. 
The positions of the more than 80 Ha emission 
sources identified here are shown in Figure 2. 

In Table 1 ordered by right ascension, we 
present identifiers for the detected Ha emis¬ 
sion sources, J2000 coordinates, spectral types, 
H—band magnitudes, V — Rc and V — Ic colors; 
J—H and H—Ks colors, and Ks magnitudes from 
2MASS; [3.6], [4.5], [5.8], and [8.0] SpitzerlRAC 
photometry from Stelzer & Scholz (2009) and ref¬ 
erences therein as well as Wide-Field Infrared Sur¬ 
vey Explorer (WISE) wl, w2, w3, and wA photom¬ 
etry. The WISE images of each source were ex¬ 


amined individually in all passbands to ensure the 
sources were valid detections, particularly in the 
confusion-limited cluster center where nebulosity 
dominates infrared emission. For the previously 
unknown emission-line sources, an IHa number 
is assigned, continuing the numbering convention 
initiated by Herbig (1998) in IC348 and subse¬ 
quently continued by Herbig & Dahm (2002) in 
IG5146, Herbig et al. (2004) in NGG1579, Dahm 
& Simon (2005) in NGG2264, Dahm (2005) in 
NGC 2362, Herbig & Dahm (2006) in L988, and 
Dahm et al. (2012) in IC 1274. The measured 
equivalent widths of the Ha emission profiles, 
W{Ha), and of Ca II A8542 are also provided 
in Table I (positive values indicate emission). 
No corrections have been applied to W{Ha) or 
W (8542) for underlying absorption structure. 

The equivalent width of Ha is a well-established 
indicator of accretion processes and chromospheric 
activity in pre-main sequence stars (for a review, 
see Bertout 1989). Traditionally the boundary 
separating classical and WTTS was placed at 
W(Ha)=10A (e.g. Herbig 1998). While no physi¬ 
cal interpretation was intended for this value, clear 
differences in the processes responsible for emis¬ 
sion have since been recognized for CTTS, i.e. ac¬ 
cretion, and WTTS, i.e. enhanced chromospheric 
activity. Various spectral type dependent criteria 
have been suggested to better distinguish accretors 
from non-accretors, e.g. Martin (1998) and White 
& Basri (2003). In Figure 13 we plot W(Ha) as a 
function of V — Ic color (top panel) and W(IIa) 
as a function of Ks — [4.5] color (bottom panel). 
If Spitzer [4.5] photometry were not available for 
a given source, Ks — w2 values are substituted. 
While no obvious correlation is present between 
W{Ha) and V — Ic color, a log-linear relationship 
is evident when comparing W(Ha) and Ks — [4.5] 
color, both disk indicators. The normal K — w2 
color for young, 5-30 Myr M5-type stars (the lat¬ 
est spectral identified in the cluster, with one ex¬ 
ception) is 0.43 mag (Pecaut & Mamajek 2013). 
Sources lying redward of this value in Figure 13 
should be considered as strong disk candidates. 

The primary sources of contamination among 
the Ha emission sources are active late-K and 
M-type field dwarfs (dMe) that exhibit enhanced 
chromospheric activity. Ha emission strengths 
among dMe stars, however, are typically weak, 
W(Ha)< 10 A (Hodgkin et al. 1995; Reid et al. 
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1995; Hawley et al. 1996), and would predom¬ 
inantly affect the statistics of the WTTS popu¬ 
lation. Other potential sources of contamination 
include chromospherically active giants, RS CVn 
binaries, cataclysmic variables, and active galax¬ 
ies. The field density of such objects is expected 
to be low, particularly within the molecular clouds 
where extinctions reach Ay ^15-20 mag. The dark 
nebulosity likely obscures background sources over 
a substantial fraction of the survey area. With 
these caveats we assume that field contamination 
among the Ha emission sources is low, but cer¬ 
tainly non-zero. 

4.3. Chandra ACIS X-ray Detected Sources 

Stelzer & Scholz (2009) obtained a 22 ks long in¬ 
tegration of NGC 7129 using ACIS onboard Chan¬ 
dra, detecting 59 X-ray sources. The majority 
(47/59) of these sources have 2MASS near-infrared 
counterparts, but prior to this survey few had op¬ 
tical photometry available. Correlating these X- 
ray detections with 2MASS and Spitzer infrared 
excess sources, Stelzer & Scholz (2009) identified 
one Class O/I source, 16 Class H sources, and 
30 Class HI candidates, leaving 12 unclassified. 
Contamination of the X-ray selected sample from 
field interlopers, particularly among the Class HI 
sources was considered, however, using X-ray lu¬ 
minosities of field dwarfs, Stelzer & Scholz (2009) 
estimate that <10“^ X-ray emitting field stars lie 
within the central core of the cluster. A contribu¬ 
tion of extragalactic X-ray sources was also consid¬ 
ered, but reduced sensitivity at large off-axis an¬ 
gles should limit the total number of extragalactic 
sources considerably. 

The X-ray detections from the Chandra ACIS 
integration are plotted as green crosses in Figure 
2. Over half, ~30, of the X-ray detections have 
counterparts in the Ha emission selected sample. 
About a half-dozen X-ray detections in the halo 
of the star forming region lacking optical counter¬ 
parts are possible extragalactic sources. Others 
appear to be associated with stars that were not 
identified as Ha emitters, but that could be clus¬ 
ter members. Clearly X-ray and Ha emission are 
tracing similar activity in these pre-main sequence 
candidates. There are a substantial number of Ha 
emission sources, however, that were not detected 
by the X-ray survey, particularly off-axis. This 
is likely the result of the reduced sensitivity and 


the relatively shallow ACIS exposure. Included in 
Table I, ordered by right ascension, are 46 of 59 X- 
ray sources identified by Stelzer & Scholz (2009), 
their J2000 coordinates, spectral types, optical {V, 
V — Rc, V — Ic) and infrared photometry from 
2MASS {J-H, H-Ks, Ks), Spitzer, and WISE. 
Of the remaining 13 sources, eight were outside of 
the fields of view of the optical photometric sur¬ 
veys and five had no optical counterparts. 

4.4. Infrared Excess Sources not Detected 
by Chandra or the Ha Emission Survey 

The Spitzer IRAC and MIPS photometry of 
Cutermuth et al. (2004) and Muzerolle et al. 
(2004) were used by Stelzer & Scholz (2009) to 
identify sources with infrared excess emission at¬ 
tributable to circumstellar disks. Stelzer & Scholz 
(2009) made their source selection using the [3.6] — 

[4.5] , [5.8] — [8.0] and the Ks — [4.5], J — H color- 
color diagrams, resulting in the identification of 
64 Class H sources and 13 Class O/I candidate 
cluster members. Of these 77 sources, 46 were 
within the Chandra field of view, but were not de¬ 
tected in X-rays. Eliminating 15 of these sources 
that were detected with Ha emission and that al¬ 
ready appear in Table 1, we list in Table 2 the 
remaining 31 pre-main sequence candidates that 
were selected on the basis of their infrared photom¬ 
etry. Contamination of these pre-main sequence 
candidates is dominated by extragalactic sources, 
predominantly star forming galaxies. Such ob¬ 
jects, however, reside in a specific region of the 

[3.6] — [4.5], [5.8] — [8.0] color-color diagram and 
can be distinguished from pre-main sequence stars 
by establishing a brightness cutoff determined by 
statistical means (Cutermuth et al. 2008; Stelzer 
& Scholz 2009). In Table 2 we present optical 
(Ic, Rc — Ic) photometry for these Spitzer excess 
sources and infrared photometry from 2MASS, 
Spitzer, and WISE. These sources are shown in 
Figure 2 as open red squares and appear to be 
preferentially positioned within the semi-circular 
arc of remnant molecular gas to north, east and 
south of the cluster core. The Class I source adja¬ 
cent to LkHa 234 is 2MASS J21430696-f660641.7, 
which is discussed in more detail in §6.2. 
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4.5. Spectroscopically Classified Sources 
Lacking Ha or X-ray Emission 

In Table 3 we present optical and infrared pho¬ 
tometry from 2MASS and WISE for sources that 
were classified by the low-dispersion spectroscopic 
survey, but that lack Ha or detected X-ray emis¬ 
sion. In §5.4, the near-infrared color-color dia¬ 
grams are used to identify several of these stars 
as possible infrared excess sources. The majority, 
however, are probable field interlopers, not associ¬ 
ated with the young cluster or the molecular cloud 
complex, but are presented here for completeness. 

5. Cluster Properties 

5.1. Reddening and Extinction 

Extinction and CO column density are in gen¬ 
eral strongly correlated, implying that dust is well- 
mixed with the molecular gas in molecular cloud 
complexes. Shown in Figure 14 is a map of vi¬ 
sual extinction {Ay) derived from integrated 
line intensity obtained at the Five College Radio 
Astronomical Observatory (FCRAO) in 1993 by 
LAH. The map is superimposed upon a near in¬ 
frared mosaic of the region with extinction peaks 
that are coincident with the compressed ridge to 
the north, east and south of the cluster core, peak¬ 
ing near Ay~20 mag. The evacuated cavity is 
clearly evident with extinction dropping sharply to 
the west and gradually to the east across the sur¬ 
face of the molecular cloud. Also shown in Figure 
14 is a map of integrated CS line intensity. The CS 
extinction map demonstrates significant enhance¬ 
ments near the molecular outflow originating by 
LkHa 234 and to the south near FIRS 2. 

The low-dispersion spectroscopy allowed spec¬ 
tral classification for ~130 stars in the cluster re¬ 
gion. Combined with the optical photometry, the 
spectral types permitted an independent determi¬ 
nation of extinction across the region. To estimate 
extinction for sources of known spectral type, we 
assume the standard ratio of total-to-selective ab¬ 
sorption, i.e. R = Ay/E{B — V)= 3.08 (He et 
al. 1995) to derive a normal reddening law given 
by Ay = 2A3E{V — Iq)- Intrinsic colors of 5- 
30 Myr old pre-main sequence stars taken from 
Pecaut & Mamajek (2013) were used to determine 
V — Ic color excesses and extinctions. The average 
extinction suffered by 50 probable cluster mem¬ 


bers taken from the Ha and X-ray selected stellar 
samples is Ay=1.8 mag. An abnormal extinction 
law induced by dust grains having sizes substan¬ 
tially larger than interstellar grains or the peculiar 
spectral energy distributions for these pre-main se¬ 
quence members could impact this adopted mean 
extinction value. 

5.2. The Color-Magnitude Diagram of the 
Cluster Population 

Shown in Figure 15 is the observed V — Ic,V 
color-magnitude diagram for all sources detected 
by the KPNO T2KA CCD imaging survey. The 
candidate members of NGC7129, i.e. the Ha 
emission and X-ray detected sources, are super¬ 
posed in the hgure. The zero age main sequence 
(ZAMS) of Siess et al. (2000) is overplotted us¬ 
ing the dwarf colors presented by Kenyon & Hart¬ 
mann (1995) and assuming a distance of 1150 pc 
(Straizys et al. 2014). The cluster sequence is well- 
dehned by the activity-selected sample and lies at 
least two magnitudes above the ZAMS, even for 
the early-type members. 

The X-ray and Ha emission populations over¬ 
lap considerably and appear to have identified 
the majority of possible cluster stars in the color- 
magnitude diagram. A collection of faint (V > 18 
mag) Ha emission sources are evident outside of 
the nominal cluster sequence lying on or just be¬ 
low the ZAMS. Most of these sources lie on the 
periphery of the molecular clouds, three to the 
west where extinction is minimal. If taken at face 
value, the isochronal ages assigned to these sources 
would be much greater than that of the cluster it¬ 
self. Such objects are generally interpreted as hav¬ 
ing edge-on disk geometries that dramatically re¬ 
duce stellar luminosity, or as background contam¬ 
inants. The roughly 10% fraction of sources lying 
substantially below the cluster pre-main sequence 
and potentially having edge-on disk geometries is 
consistent with having a random distribution of 
disk orientations. 

In Figure 16 we show the extinction-corrected 
(V — Ic)o, Vo color-magnitude diagram for the Ha 
emission stars and X-ray sources. Stars of known 
spectral type have been corrected individually for 
reddening using the intrinsic colors of 5-30 Myr 
pre-main sequence stars from Pecaut & Mamajek 
(2013). All other sources have been dereddened 
using the mean extinction value derived for the 
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Ha and X-ray samples. 

5.3. The Ages and Masses of the Pre-Main 
Sequence Population 

The age of the cluster emerging from NGC 7129 
is undoubtedly young given the presence of one 
and possibly two Herbig Be stars (LkHa 234, SVS 
13) as well as the classical Be star BD-|-65°1637. 
As will be discussed in §6, the evolutionary 
timescale of the photodissociation region envelop¬ 
ing BD-|-65°1638 may be even younger, <10^ yr. 
Cluster ages are generally inferred by fitting their 
stellar populations in color-magnitude diagrams 
with pre-main sequence isochrones. 

The grid of Siess et al. (2000) models was used 
to estimate the ages and masses for the optically 
detected Ha emitters and X-ray sources having 
established spectral types. Uncertainties involved 
in the use of pre-main-sequence models fall into 
two broad categories; the physics used in model¬ 
ing stellar evolution from the birthline to the zero- 
age main sequence, and the transformation be¬ 
tween theoretical and observational planes. Mod¬ 
els of pre-main sequence evolution treat convec¬ 
tion, opacity, radiative transfer, rotation, and ac¬ 
cretion uniquely, leading to variations in predicted 
evolutionary paths for a given stellar mass. Initial 
conditions establishing the birthline are also not 
well understood. For an in depth discussion of is¬ 
sues related to the ages of young stars, the reader 
is referred to Soderblom et al. (2014). 

Transforming between theoretical and obser¬ 
vational planes is typically achieved by fitting 
main-sequence colors and bolometric corrections 
as a function of effective temperature, a problem¬ 
atic assumption given the lower surface gravities, 
cooler temperatures, enhanced chromospheric ac¬ 
tivity and accretion processes generally associated 
with pre-main-sequence stars. The Siess et al. 
(2000) models adopt the dwarf colors of Kenyon 
& Hartmann (1995), which themselves are derived 
from Bessell & Brett (1988). 

With these caveats stated, the median age for 
the 50 cluster members having established spec¬ 
tral types determined from the models of Siess 
et al. (2000) is ~1.8 Myr. The majority of 
sources in Figure 16 fall between the 1 and 3 
Myr isochrones, implying relative coevality, within 
the tolerance of expected errors that are summa¬ 


rized by Soderblom et al. (2014). The spread in 
extinction-corrected, visual magnitude for a given 
color or effective temperature varies between one 
and two mag for U — Ic >l-5- A literal reading 
of the implied distribution of ages suggests a full 
range from 0.1 Myr to ~35 Myr; however some of 
the more advanced ages are potentially attributed 
to neutral extinction induced by edge-on disk ge¬ 
ometries. These ages are sensitive to the assumed 
distance of the cluster, intrinsic variability, errors 
in extinction, as well as unresolved binaries. The 
extreme case of having companions of equal mass 
would elevate stars in the color-magnitude dia¬ 
gram by 0.75 mag, resulting in the assignment of 
younger ages. 

Recently, revisions to age estimates of young 
clusters and associations have advanced published 
ages by a factor of two or more. Pecaut et al. 
(2012) fit isochrones to the F and G-type stars in 
the color-magnitude diagram of the Upper Scor- 
pius OB association, suggesting a revision to its 
age from ^5 Myr, which characterizes the K and 
M-type stellar population, to 11 Myr if effects of 
the stellar birthline in the HR diagram (resulting 
from the finite radii of collapsing protostars) can 
be neglected. Bell et al. (2013) fit cluster popu¬ 
lations with new semi-empirical isochrones which 
adopt empirically derived color-effective tempera¬ 
ture relationships and bolometric corrections for 
pre-main sequence stars. The resulting ages of 
several young clusters and star forming regions in¬ 
cluding IC5146, IC348, and NGG2362 were sub¬ 
sequently adjusted by Bell et al. (2013) from <1 
Myr to ~2 Myr, from ~1.3 Myr to ~6 Myr, and 
from ~5 Myr to ~12 Myr, respectively. Appli¬ 
cation of such techniques to the population of 
NGG 7129 would similarly advance its age from 
the estimate derived here if these techniques and 
assumptions are robust to further tests. 

The inferred masses of the cluster members 
range from ~0.2 Mq to 5.6 Mq (the most mas¬ 
sive being BD-|-65°1638), implying under standard 
initial mass function assumptions that a substan¬ 
tial fraction of low-mass stars and brown dwarfs 
remain unaccounted for in the activity-selected 
sample. Assuming a photometric completeness 
limit of V~21 within the central cluster region im¬ 
plies that the survey is complete to about ~0.3 
Mq. The Keck LRIS Rc and Ic—band imag¬ 
ing is the deepest of the three photometric sur- 
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veys presented here, extending to /c~22 within 
the core of NGC 7129 where extinction is signif¬ 
icantly reduced. Shown in Figure 17 is the ob¬ 
served Rc — Ic) Ic color-magnitude diagram for 
the Keck LRIS photometry, uncorrected for red¬ 
dening. The Siess et al. (2000) models do not 
extend below 0.1 M©, but it is clear from the fig¬ 
ure that about two dozen sources are present that 
would lie along the 2 Myr isochrone if extrapo¬ 
lated to lower masses. The evolutionary models of 
Baraffe et al. (1998) extend to much lower masses, 
well below the substellar mass limit and into the 
planetary mass regime. Assuming a mean extinc¬ 
tion of Av ~1.8 mag or Ai ~1.0 mag, the Baraffe 
et al. (1998) models predict that the sub-stellar 
mass limit for NGC 7129 lies near Mj = 8.53 or 
/c=19.83 mag in Figure 17. Here there are a 
handful of sources in the photometric cluster se¬ 
quence that are possible very low mass stars or 
brown dwarf candidates, but a deeper photomet¬ 
ric survey as well as spectroscopy, particularly in 
the near-infrared, are needed for confirmation. 

5.4. The Infrared Color-Color Diagram 
for the Cluster Population 

NGC 7129 was observed from 3.6 to 24 /xm 
with Spitzer using IRAC and MIPS as part of 
the Spitzer Young Stellar Cluster Survey (AOR: 
3655168 and 3663616). These observations are 
presented by Megeath et al. (2004) as well as 
Gutermuth et al. (2004, 2009). Stelzer & Scholz 
(2009) adopted the photometry of Gutermuth et 
al. (2004, 2009) to derive a disk fraction for the 
Chandra X-ray and infrared-selected sources of 
~51±14%. In a mass-limited, lightly extincted 
sub-sample (Ay < 5 mag) of pre-main sequence 
stars, Stelzer & Scholz (2009) derive a disk frac¬ 
tion of ~33ti9%. Although considered the least- 
biased sub-sample, it is also the least complete in 
terms of the number of included cluster members. 

Here we combine 2MASS, Spitzer IRAC, and 
WISE photometry for the activity-selected sources 
in NGC 7129 to examine infrared excesses for can¬ 
didate cluster members. Shown in Figure 18 are 
the H — Ks, J — H (left panel) and the ATg — [4.5], 
J — H (right panel) color-color diagrams for the 
Ha emission sources. X-ray sources, infrared ex¬ 
cess sources not identified with X-ray or Ha emis¬ 
sion, and the spectroscopically classified sources 
lacking either Ha or X-ray emission. If IRAC [4.5] 


channel photometry were not available, WISE w2 
photometry was substituted. The main sequence 
colors of Pecaut & Mamajek (2013) are overplot¬ 
ted and the approximate reddening boundaries for 
dwarfs are shown with slopes derived using the ex¬ 
tinction coefficients for diffuse interstellar clouds 
taken from Martin & Whittet (1990). The 2MASS 
H — Ks, J — H color-color digram distinguishes 
the most prominent infrared excess sources that 
lie to the right of the reddening boundary for nor¬ 
mal dwarfs. These excesses arise from hot >1000 
K dust lying within the innermost (0.1-1 AU, de¬ 
pending upon stellar luminosity) regions of the cir- 
cumstellar disks. The Ks — [4.5], J — H color- 
color diagram is much more effective at isolating 
disk-bearing sources given that the IRAC [4.5] and 
WISE w2 photometry is sensitive to cooler dust 
temperatures. 

Assuming the reddening boundary for late-type 
dwarfs in the Ks — [4.5], J — H color-color diagram 
to represent the demarcating limit for disk-bearing 
sources, we find the disk fraction for the activity- 
selected sample having [4.5] or w2 photometry 
available to be 57±9%, consistent with the disk 
fractions derived by Gutermuth et al. (2004), 
54±I4%, and that of Stelzer & Scholz (2009), 
~5I±14%. This should be considered as a lower 
limit for the sample examined here given that stars 
having weak excesses and lying within the color 
boundaries of normal stars have been excluded 
from the disk population. The estimated disk frac¬ 
tion in NGC 7129 is similar to those derived for 
the comparably aged star forming regions IC 348, 
50±6% (Lada et al. 2006) and Chamaeleon I, 
50±6% (Luhman et al. 2008), but substantially 
larger than the disk fractions of the more evolved 
Upper Scorpius OB association, 19% (Carpenter 
et al. 2006) and NGC 2362, 19% (Dahm & Hillen- 
brand 2007). 

Eigure 18 suggests that Ha emission is a better 
discriminant for identifying disk-bearing sources 
than X-ray emission. Longer exposure times with 
ACIS onboard Chandra would certainly have iden¬ 
tified additional cluster members, but the effec¬ 
tiveness of simple ground-based, slitless grism or 
multi-object spectroscopic techniques at isolating 
active stars cannot be ignored. It should also be 
noted that a handful of spectroscopically classified 
sources lacking either Ha or X-ray emission appear 
to have substantial infrared excess and are there- 
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fore candidate cluster members. These sources are 
identified as having infrared excess in the com¬ 
ments column of Table 3. The majority of the 
classified sources lacking activity indicators, how¬ 
ever, appear to cluster along the upper extinction 
boundary for normal dwarfs where reddened gi¬ 
ants might be expected to lie, supporting their 
non-membership status in NGC 7129. In sum¬ 
mary, to conduct a population census of young 
star forming regions, Ha emission. X-ray emission 
and infrared excess are needed to ensure complete¬ 
ness. 

6. Discussion 

6.1. Star Formation in NGC 7129 and the 
105?4-|-09?9 Molecular Cloud Com¬ 
plex 

Star formation is at an advanced stage within 
NGC 7129 with the recently formed massive stars 
BD-|-65°1638 and BD-|-65°1637 dissociating and 
dispersing molecular material in the center of the 
molecular cloud. The exact sequence of massive 
star formation has been debated in the literature 
with Bechis et al. (1978) and Miskolczi et al. 
(2001) suggesting that BD-|-65°1638 collapsed ini¬ 
tially and that its ultraviolet flux alone is largely 
responsible for clearing the elongated cavity. 

Matthews et al. (2003), however, argue that 
winds from BD-|-65°1637 excavated the cavity 
given that BD-|-65°1637 is the more luminous of 
the pair and that in projection, is more centrally 
located within the cleared area. Their moderate- 
resolution (1') 21 cm H I survey of NGC 7129 
revealed a ~30' diameter ring of emission appar¬ 
ently associated with the surface of the molecular 
cloud. H 2 on the cloud surface has likely been 
dissociated by interstellar ultraviolet radiation, 
producing the observed ring. Within this ring, 
Matthews et al. (2003) note a bright knot of H I 
emission centered on BD-|-65°1638 that extends 
well beyond the confines of the molecular ridge. 
They suggest that the H I knot and BD-f65°1638 
are seen in projection against the near surface of 
the molecular cloud, and are not lying within it. 
This would account for the observed radial veloc¬ 
ity difference of ~2 km s“^ between the CO and 
the H I emission. Matthews et al. (2003) postu¬ 
late that BD-|-65°1638 is the younger of the two 
stars having only formed within the last 1500-10^ 


yrs on the cloud periphery. This would also im¬ 
ply that BD-|-65°1638 is among just a handful of 
“dissociating stars” known, massive stars caught 
in the process of emerging from their natal shells. 

The HIRES spectra of BD-t65°I637 and BD-t65°I638 
clearly demonstrate that BD-|-65°1637 is the more 
active of the pair. Strong Ha emission is apparent 
in the former, presumably formed by recombina¬ 
tion within the surrounding gaseous disk. Stellar 
winds are responsible for the blue-shifted forbid¬ 
den emission present in the spectrum. In contrast, 
the spectrum of BD-|-65°1638 is unremarkable 
with little evidence to suggest the extreme youth 
required of it to be a dissociating star. In the 
infrared, strong excess emission characteristic of 
circumstellar material does not become apparent 
in the spectral energy distribution of BD-|-65°1638 
until ~8 ^m. The star falls on the locus of the 
main sequence in the near infrared Ks — [4.5], 

J — H color-color diagram shown in Figure 18. 
BD-|-65°1637 lies well to the right of the sequence 
suggestive of infrared excess. BDH-65°1637 and 
BD-f65°1638 lie near each other on the extinction- 
corrected, color-magnitude diagram shown in Fig¬ 
ure 16 and suffer comparable levels of optical ex¬ 
tinction, Ay =1.7 and 2.1 mag, respectively. 

The preponderance of spectroscopic evidence 
suggests that BD-|-65°1638 is the more quiescent 
and therefore more evolved star of the pair. The 
origin of the bright knot of H I emission is un¬ 
resolved, but may represent an expanding shell 
of dissociated molecular gas driven by the UV 
flux from BD-t65°1637 and BD-f65°1638. The 
presence of an extensive photodissociation region 
(PDR) inferred from an arc of H 2 emission some 
170" long or 0.95 pc, was noted by Schultz et al. 

(1997) in their narrow-band, near-infrared imag¬ 
ing survey of NGG7129. This PDR appears to 
outline the molecular ridge and is thought to be 
illuminated by UV emission from BD-|-65°1638. 
Regardless of formation order, these early-type 
stars have disrupted the star formation process 
within LI 181, effectively dispersing molecular ma¬ 
terial within the cavity and possibly induced star 
formation throughout the molecular cloud. 

Triggered star formation has been suggested for 
the central cluster region (e.g. Miskolczi et al. 

2001), particularly along the dense molecular ridge 
where LkHa 234 and 2MASS J21430696-h660641.7 
are located. Maps of GO column density (e.g. 
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Bechis et al. 1978; Miskolczi et al. 2001; Figure 
14a), reveal another dense core ~1 pc south of 
BD+65°1638 where the intermediate-mass Class 
0 protostar FIRS 2 is located. The mass of this 
source has been estimated by Fuente et al. (2014) 
to be ~2-8 M 0 . Whether triggered by com¬ 
pression from the formation of BD-|-65°1637 and 
BD-|-65°1638 or by some other event is unknown. 
Surprisingly few Ha emission sources or embedded 
Class 0 and I sources are found in this region, pos¬ 
sibly implying that star formation is in its earliest 
phase here. 

The activity-selected sample of candidate mem¬ 
bers is nearly evenly divided between those con¬ 
centrated within the cluster core, defined as lying 
within 2' (~0.7 pc) of BD-|-65°1638, and a dis¬ 
persed population spread throughout the molecu¬ 
lar cloud complex. The majority of this dispersed 
population lies within the northeast quadrant of 
L1181, extending from V350 Cep and SVS 10 to 
south of GGD 35 (HH 235). This broad arc spans 
more than three parsecs in projection and is char¬ 
acterized by the presence of HH objects, shocked 
molecular hydrogen emission, and molecular out¬ 
flows. Over three dozen Ha emission sources, X- 
ray sources and embedded protostars have been 
identified, suggesting that star formation is well 
underway in this region of the molecular cloud. 
These isolated sources have likely formed indepen¬ 
dently of the central young cluster. 

To the west, several Ha emitters are located 
around the molecular cloud core GGD 32, with HH 
103 marking the terminus of star formation activ¬ 
ity. No embedded protostars are found here and 
the detected Ha sources are predominantly knots 
of HH objects. The paucity of activity-selected 
sources further west suggests that few have drifted 
far from their formation sites, arguing for recent 
and rapid star formation in L1181. 

6.2. The Origin of the Molecular Outflow 
Near LkHa 234 

LkHa 234 lies at the apex of the one-parsec 
scale, central cavity, the boundaries of which are 
demarcated on optical and infrared continuum im¬ 
ages by bright arcs of nebulous emission and an 
increased density in stars, both cluster members 
and background sources. Northeast of NGG 7129, 
bright 4.5 /im emission is evident in the Spitzer 
IRAG 4-color image of the region presented by 


Gutermuth et al. (2004). This shocked GO emis¬ 
sion is possibly produced by the impact of super¬ 
sonic jets originating near LkHa 234. LkHa 234 is 
also a radio continuum source (Bertout & Thum 
1982; Snell & Bally 1986) and was thought to 
power an optical jet (Ray et al. 1990). The calcu¬ 
lated mass loss rate for the optical outflow is ~ 2 - 
3x10“® Mq yr“^ and implies a momentum flux 
which is two orders of magnitude less than that 
of the molecular outflow. The optical jet points 
along the symmetry axis of the molecular cavity, 
with a number of faint stars aligned along the same 
axis. Because of its high collimation, however, the 
jet is unlikely to have formed the parsec-scale cav¬ 
ity. Ray et al. (1990) report that the HH objects 
present throughout the region are not associated 
with LkHa234 itself based upon a proper motion 
analysis. Shocked H 2 emission was reported by 
Wilking et al. (1990), Schultz (1989), and Schultz 
et al. (1997) south of LkHa234. This emission 
which could be related to the optical jet and to 
the GO molecular outflow is located on the oppo¬ 
site side of the star. 

To investigate the source of the optical jet and 
the molecular outflow near LkHa 234, NIRG2 on 
Keck H imaged the region using its wide camera in 
J, H, and K'. Shown in Figure 19 is a median com¬ 
bined, ~30"x30" K' image of the region having 
an effective integration time of 100 s. The arc of 
nebulosity evident in the NIRG2 image represents 
the apex of the ^170" long rim of H 2 emission dis¬ 
cussed by Schultz et al. (1997). LkHa 234 appears 
to be involved with this nebulosity, blowing back 
the gas and dust to form an illuminated wedge- 
shaped rim. The AO imaging reveals at least a 
dozen sources including a well-resolved companion 
of LkHa 234 having a position angle (PA) of 97° 
and a separation of 1788. This source was noted 
previously by Perrin (2006) in near-infrared polari- 
metric observations of the region. About ~12'.'5 
south of LkHa 234 lies the embedded Class I pro¬ 
tostar 2MASS J21430696-I-660641.7, which only 
becomes prominent in the NIRC2 imaging at K'. 

The NIRC2 imaging also reveals two point 
sources along the southern arm of the nebulous 
arc, including a resolved visual binary with a sep¬ 
aration of ^0733 and AmKp~2.1 mag. A lumi¬ 
nous knot of nebular emission is apparent 277 dis¬ 
tant from LkHa 234, PA 280°. Weintraub et al. 
(1994) reported an embedded companion lying 3" 
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northwest of LkHa 234 (PA 290-340°), their po¬ 
larization source PS 1, and suggested that this 
source is responsible for driving the bipolar molec¬ 
ular outflow associated with the region. PS 1 
and their source IRS 5 (PA=330°, 2"5 separa¬ 
tion from LkHa234) roughly correspond with this 
bright knot of nebulous emission. The axis of the 
primary molecular outflow in the region is approx¬ 
imately 60°/240°, which corresponds with a sub¬ 
stantial number of HH objects identified in the 
area (McGroarty et al. 2004). 

Perrin (2006) describes two mid-infrared sources 
~3-4" northwest of LkHa 234, one triangular in 
shape and the other more point source-like. These 
sources are not detected in the NIRC2 imaging, 
but would be roughly coincident with PS 1 and 
the knot of nebulous emission. Reproduced in Fig¬ 
ure 20 is an archival Spitzer IKPiC [5.8] post basic 
calibrated data (BCD) image of the region, cen¬ 
tered near LkHa 234. The PSF core of LkHa 234 
is saturated in the image as is the Class I source 
2MASS J2I430696-I-660641.7. A bright infrared 
source corresponding to the mid-infrared sources 
of Perrin (2006) is evident as is a third point 
source ~10" northwest of LkHa 234. These in¬ 
frared sources appear to be roughly aligned with 
the bow-shocked nebular emission enshrouding 
SVS 13 further to the west. This bow-shocked 
emission either results from winds emanating from 
SVS 13 impacting surrounding molecular material 
or possibly from outflowing molecular gas from 
these embedded infrared sources colliding with 
winds from the early-type star. 

Between LkHa 234 and 2MASS J21430696-b660641 
lies a third bright infrared point source in the 
IRAC [5.8] image that lacks counterpart in the 
NIRC2 K' image. This source is readily appar¬ 
ent in Spitzer IRAC [8.0] imaging and is likely 
another embedded Class 0 or I protostar. Assum¬ 
ing all NIRC2 and IRAC sources in this region 
to be cluster members, the stellar surface density 
in the immediate vicinity of LkHa 234 is several 
hundred (>400) stars pc“^. Civen the number 
of embedded sources emerging here, this region 
at the apex of the parsec-scale cavity must be 
extremely young. There are no other comparable 
examples of clustered, luminous protostars evident 
in the region. 


7. Summary 

We have obtained deep, optical {V,Rc,Ic) 
photometry for ~2500 sources in a field roughly 
centered on the reflection nebula NCC 7129 in 
Cepheus. Over 80 Ha emission sources have been 
identified in the region, the majority of which 
are presumably members of a T Tauri star pop¬ 
ulation emerging from the associated molecular 
cloud. Combined with 59 X-ray sources detected 
in a shallow, 22 ks Chandra ACIS observation by 
Stelzer & Scholz (2009), the Ha emission sources 
form a relatively narrow pre-main sequence in the 
V — Ic, V color-magnitude diagram. Inspection of 
the color-magnitude diagram for all sources sug¬ 
gests that these two stellar activity indicators have 
identified the majority of the optically-detectable 
pre-main sequence population, down to a limiting 
magnitude of y~20. 

Spectral types for over 130 sources in the region 
have been determined from low-dispersion, red op¬ 
tical spectroscopy. The mean extinction suffered 
by a sample of candidate cluster members is about 
Ay ~1.8 mag. Using the pre-main sequence evo¬ 
lutionary models of Siess et al. (2000), we derive 
a median age of ~1.8 Myr for the Ha and X-ray 
selected sources having established spectral types, 
although a substantial age dispersion is present 
(~0.1-35 Myr). 

Using near and mid-infrared photometry from 
2MASS, Spitzer and WISE, we estimate the disk 
fraction of the activity- selected sample of clus¬ 
ter members to be ~57±9%. A small fraction of 
spectroscopically classified stars lacking either Ha 
or X-ray emission that exhibit infrared excess are 
present. We find a strong log-linear relationship 
between W(Ha) and Ks — [4.5] color, as might be 
expected from these two inner disk indicators. 

Merging the activity-selected sample of 94 
sources with 31 infrared selected pre-main se¬ 
quence candidates from Stelzer & Scholz (2009) 
not detected with X-ray or Ha emission, and eight 
infrared excess sources from the spectroscopically 
classified sample of stars lacking X-ray or Ha emis¬ 
sion, the cluster population of NCC 7129 has been 
increased from ~90 members to more than 130. 

High dispersion optical spectroscopy of BD-|-65°1638 
and BD-|-65°1637 obtained with HIRES on Keck 
I are presented which clearly demonstrate that 
BD+65°1637 is the more active of the two, a 
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probable classical Be star. The spectrum of 
BD+65°1638 is devoid of emission, suggesting 
the star to be remarkably quiescent if indeed a 
photodissociating star, just emerging from its pri¬ 
mordial envelope of molecular gas. 

HIRES spectroscopy of LkHa 234 is suggestive 
of active accretion with numerous metallic and for¬ 
bidden transitions in emission as well as Ha and 
H/3. HIRES spectra obtained in 1999 and 2004 
show significantly different Ha emission profiles 
including the presence of a deep absorption fea¬ 
ture in the P Cygni-like structure in the earlier 
spectrum. The high dispersion spectrum of SVS 
13 reveals weak emission reversal in the core of 
Ha as well as Ee I AA4988, 4994 and [O I] A6300 
emission. Three of four early-type stars within the 
cluster core are confirmed as Be stars. 

Adaptive optics imaging of LkHa 234 using 
NIRC2 on Keck H reveals a number of embedded 
sources in the region including the Class I proto¬ 
star 2MASS J21430696-b660641.7. A bright, well- 
resolved companion of LkHa 234 is also detected 
having a PA of 97° and a separation of 1'.'88. Of 
note several mid-infrared sources that are appar¬ 
ent in Spitzer IRAC and MIPS imaging of the re¬ 
gion are not detected in the Keck AO K' imaging. 
Assuming all NIRC2 and IRAC sources to be asso¬ 
ciated with the cluster, the stellar surface density 
near LkHa 234 is several hundred stars pc“^. 

Our interpretation of the structure of the rem¬ 
nant molecular clouds and the distribution of 
young stars in the region is that BD-|-65°1638 
is primarily responsible for evacuating the parsec- 
scale cavity in LII8I. LkHa 234 and a small 
number of deeply embedded protostars are form¬ 
ing within the ridge of compressed gas, possi¬ 
bly triggered by the collapse of BD-|-65°I638 and 
BD-|-65°1637. Star formation within the compact 
cluster has been rapid and coeval given the rel¬ 
atively narrow pre-main sequence evident in the 
color-magnitude diagram. The formation of other 
intermediate-mass stars to the north and south 
of the cluster core including SVS 10 and FIRS 
2, will likely further disrupt the molecular cloud, 
effectively ending the era of star formation in the 
region. 

This paper is dedicated to the memory of 
George H. Herbig who passed away on 12 Oc¬ 
tober 2013 at the age of 93 after a career spanning 


more than 70 years. He is considered by many to 
be the father of observational star formation, and 
as a testament to the impact of his contributions 
upon the field, classes of pre-main sequence ob¬ 
jects now bear his name. An intellectual giant, his 
extraordinary physical insight, natural curiosity, 
and uncanny ability to ask the right questions led 
to tremendous advances in our understanding of 
the star formation process. There are few papers 
today on the subject of star formation that can¬ 
not trace their origins back to work initiated by 
George during his time at Lick Observatory, the 
University of California at Santa Cruz, and at the 
Institute for Astronomy, University of Hawaii. 
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Fig. 1.— yi^c’.^c—band three-color image of 

NGC 7129 and the L1181 molecular cloud ob¬ 
tained by G. H. Herbig in 1999 October (orien¬ 
tation: north up, east left). This 7f5x7'5 field is 
centered on the grouping of early-type stars within 
the reflection nebula that includes BD-|-65°1638 
(south), BD-|-65°1637 (west), LkHa 234 (east) and 
SVS 13 (north). LkHa 234 lies within the dense 
ridge of molecular gas that extends several ar- 
cminutes and outlines the evacuated region of the 
molecular cloud. To the northeast, just over half 
the distance from the cluster center to the corner, 
is the ghostly outline of HH 105, and to the south¬ 
west, near the edge of the field of view, HH 103 
which is likely associated with FIRS 2 or V392 Cep 
(RNO 138). FIRS 2, undetected at optical wave¬ 
lengths, lies near the center of the bottom edge of 21 
the field of view, just south of V392 Cep. 
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Fig. 2.— A 20'x20' field of the red DSS image 
of NGC 7129. The 150 square arcminute area sur¬ 
veyed with the WFGS and the region imaged with 
LRIS are outlined. The KPNO T2KA imaging 
data encompasses the entire field of view. Filled 
red circles mark the positions of Ha emission 
sources, green crosses represent X-ray detected 
sources from the Chandra AGIS observation, and 
open red squares show the positions of Glass 0/1/11 









V 


hJ M ^ ^ ^ 

NJ o CO cn 4^ ho o 



CJ)-l^h000001>NJ 



Fig. 3.— (left panel) The observed V — Iq, V 
color-magnitude diagram for all ~2500 sources de¬ 
tected iu the KPNO 1993 CCD imaging survey of 
a 23'x23' region, (right panel) The Rc — Ic, Rc 
color-magnitude diagram for all sources detected 
in the Keck LRIS 8' x 6' imaging survey of the core 
of NGC 7129. No corrections for interstellar red¬ 
dening have been applied to the photometry. 
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Fig. 4.— Comparison of the KPNO photometry 
from 1993 June with that from the UH 2.2 m from 
1999 October. The slight offset present in the V — 
Ic colors arises from extended red transmission in 
the UH 2.2 m Ic filter. 
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Fig. 5.— Comparison of the KPNO photometry 
from 1993 June with that from the Keck LRIS 
imaging survey of 1999 June. Photometric cali¬ 
bration of the LRIS data was achieved by apply¬ 
ing airmass and zero point corrections only, i.e. no 
color terms were applied, which likely explains the 
color offset. 
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Fig. 6.— Examples of low-dispersion spectra of 
late-type cluster members obtained using HYDRA 
on WIYN. The sources are ordered by spectral 
type, classified using the temperature sensitive 
indices measuring the TiO bandhead absorption 
strengths at A7140 and A8465 (light gray) relative 
to nearby continuum levels at 7035 A and 8415 A 
(dark gray). 
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Fig. 7.— Sections of the HIRES spectrum of the 
B3e star BD+65°1637 centered near Ha (left) and 
H/3 (right). Strong Ha emission {W = —25 A) is 
apparent with wings extending beyond ±550 km 
s“^. There is no indication of an underlying pho- 
tospheric absorption profile. The emission line is 
double-peaked, and the radial velocities of the red 
and blue peaks are annotated in the figure. H/3 re¬ 
veals a similar double-peaked emission structure, 
but is substantially weaker {W = —2.4 A) and re¬ 
veals the wings of a broad photospheric absorption 
prohle. 
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Fig. 8.— Sections of the HIRES spectrum of the 
B3 star BD+65°1638 centered near Ha (left) and 
H/3 (right). The spectrum reveals weak evidence 
for possible emission reversal within the core of 
the Ha absorption feature with a radial veloc¬ 
ity of 25.8 km s“^, consistent with that of 
the molecular cloud. Otherwise the spectrum of 
BD-|-65°1638 is quite unremarkable in contrast to 
those of the classical Be star BD-|-65°1637 and the 
Herbig Be star LkHa 234. 
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Fig. 9.— Sections of the HIRES spectra of LkHa 
234 obtained in 1999 December and 2004 June cen¬ 
tered near Ha showing substantial differences in 
structure between the two epochs. The 1999 spec¬ 
trum reveals a deep absorption feature in the P 
Cygni-like profile that is not evident in the later 
spectrum. A second narrow absorption core is ev¬ 
ident near line center, just blueward of the sharp 
emission peak. The emission peak and the stepped 
redward slopes are remarkably similar in both ob¬ 
servations. 
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Fig. 10.— Sections of orders from the 2004 HIRES 
spectrum of LkHa 234 centered near H/3 (left 
panel), Fe II A5018 (center panel) and the Na I 
D lines (right panel). Radial velocities of various 
features in the spectra are annotated for reference. 
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Fig. 11.— Sections of orders from the 2010 HIRES 
spectrum of SVS 13 centered near Ha (left panel) 
and H/3 (right panel). The heliocentric radial ve¬ 
locity of the emission core centered within the 
broad Ha absorption profile is annotated. SVS 
13 is possibly a Herbig Be star, but shows sub¬ 
stantially less emission than either LkHa 234 or 
the classical Be star BD-|-65°1637. 
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Fig. 12.— Sections of the HIRES spectrum of the 
enigmatic M2-type star V350 Cep centered near 
the Na I D lines (left panel), the numerous metallic 
emission lines near A6220 (center panel), and the 
P Cygni-like profile of Ha (right panel). Radial 
velocities of various features are indicated in the 
panels for reference. 
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Fig. 13.— (top panel) V — Iq color plotted against 
W(Ha) for the Ha emission sources listed in Ta¬ 
ble 1. For reference the colors of 5-30 Myr stars 
are plotted along the abscissa. No obvious corre¬ 
lation is present, but the colors shown here have 
not been corrected for extinction, (bottom panel) 
Ks — [4.5] color plotted as a function of W(Ha) 
for the Ha emission sources identified here. The 
reddest color expected from an unreddened stel¬ 
lar photosphere of M5 spectral type is indicated 
by the broken vertical line. A log-linear relation¬ 
ship is apparent between these two accretion disk 
parameters. 
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Fig. 14.— (a) A map of visual extinction (Ay) de¬ 
rived from integrated line intensity obtained 

at the Five College Radio Astronomical Observa¬ 
tory (FCRAO) in 1993 by LAH. The contour map 
is overlaid upon a near infrared mosaic image of 
NGC 7129 obtained on the KPNO 50-inch tele¬ 
scope using the SQIID camera. Extinction con¬ 
tours are plotted with intervals of 5 mag over the 
range from 2.5 to 60 mag. (6) A map of integrated 
CS line intensity from data obtained at FCRAO 
superposed upon the same near infrared mosaic 
image of NGC 7129. Extinction contours are plot¬ 
ted with intervals of Ay ~0.25 mag over the range 
from 1.0 to 2.5 mag. 
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Fig. 15.— The observed V — Ic,V color- 
magnitude diagram of NGC7129 with the Ha 
emission sources shown as solid red circles, X-ray 
sources as green crosses, and other sources from 
the KPNO T2KA survey as gray points. The 
zero age main sequence (ZAMS) of Siess et al. 
(2000) is overplotted assuming the dwarf colors 
presented by Kenyon & Hartmann (1995), derived 
from those of Bessell & Brett (1988) and a distance 
of 1150 pc (Straizys et al. 2014). 


35 




-101 234 


Fig. 16.— The extinction-corrected {V — Ic)o, 
Vo color-magnitude diagram for the Ha emission 
stars and X-ray sources with available optical pho¬ 
tometry. Stars of known spectral type have been 
corrected individually for reddening using the in¬ 
trinsic colors of 5-30 Myr pre-main sequence stars 
from Pecaut & Mamajek (2013). Sources with¬ 
out spectral type information are plotted using the 
mean extinction derived for the cluster Ay = 1.79 
mag. Superposed are the 0.5, 1.0, 3.0, and 10 
Myr isochrones and the 0.3, 0.5, and 1.0 Mq evo¬ 
lutionary tracks of the solar metallicity models of 
Siess et al. (2000), assuming a distance of 1150 pc 
(Straizys et al. 2014). 
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Fig. 17.— The observed Rc — Ic, Ic color- 
magnitude diagram of NGC 7129 constructed from 
the Keck LRIS photometry, uncorrected for red¬ 
dening. Photometric errors are shown for all 
sources. The Siess et al. (2000) isochrones for 
2 Myr, the approximate cluster age, and 100 Myr 
and the evolutionary track for a 0.1 Mq star are 
overplotted. The shaded region represents the 
area where very low mass cluster members are ex¬ 
pected to lie, including brown dwarf candidates. 
The Baraffe et al. (1998) models predict that 
the sub-stellar mass limit lies near /c=19.83 mag, 
where a handful of sources are evident in the clus¬ 
ter sequence. Deeper photometric surveys as well 
as spectroscopy are needed for confirmation. 
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Fig. 18.— The H — Ks, J — H (left panel) and 
the Ks — [4.5], J — H (right panel) color-color 
diagrams for the Ha emission sources (red cir¬ 
cles), X-ray sources (green crosses), infrared ex¬ 
cess sources lacking X-ray emission or Ha emis¬ 
sion (red squares) and spectroscopically classified 
stars lacking either Ha emission or X-ray emission 
(gray circles). The main sequence colors of Pecaut 
& Mamajek (2013) are overplotted in both pan¬ 
els as dashed lines and the approximate redden¬ 
ing boundaries for dwarfs are shown with slopes 
derived using extinction data for diffuse interstel¬ 
lar clouds from Martin & Whittet (1990). As¬ 
suming the extinction boundary for dwarfs in the 
Ks — [4.5], J — H color-color diagram to repre¬ 
sent the demarcation line for disk-bearing sources, 
we find the disk fraction for all activity-selected 
sources (i.e. X-ray, Ha emission) included here to 
be ~57±9%. 
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Fig. 19.— A median combined, ^30"x30" K' 
image centered near LkHa 234 obtained using 
NIRC2 and the Keck II NGS AO system. The 
arc of nebulosity evident in the NIRC2 image rep¬ 
resents the apex of the ^170" long rim of H 2 
emission discussed by Schultz et al. (1997). A 
well-resolved companion of LkHa 234 having a 
PA of 97° and a separation of 1"88 is evident 
as is the embedded young Class I source 2MASS 39 
J21430696-I-660641.7, some 12'.'5 to the south. 
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Fig. 20.— Spitzer IRAC [5.8] post BCD image 
centered near LkHa 234. For reference, the NIRC2 
field of view is outlined in white. LkHa 234 and 
2MASS J21430696+660641.7 are saturated in the 
image, resulting in black cores for these sources. 
Three luminous infrared sources are evident that 
were not detected in the NIRC2 imaging of the 
region. It is possible that one of these sources is 
responsible for the molecular outflow and for the 40 
bow-shocked emission enveloping SVS 13 to the 
west. 



Table 1 

Optical and Infrared Photometry for Ha and X-ray Emission Sources in NGC7129 


Identifier®' 

(J2000) 

S 

(J2000) 

SpT 

yb 

V - Be*’ 

V - /c*’ 

J - 

H - 

^3° 

[3.6]'* 

N.s]-* 

[5.8]'* 

[8.0]'* 


Tn2® 

it;3® 

it)4® 

W(H<1* 

(A) 

W(85428 ) 

(A) 

Other Identifiers 

Notes^ 

MMN 1 

21 

42 

23.08 

+66 

06 

04.4 

M3V 

20.47 

1.56 

3.11 

1.01 

0.66 

13.39 

12.51 

12.14 

11.78 

11.02 

12.39 

12.02 



32.9 

0.3 

S3-U840 


IHa 763 

21 

42 

26.60 

+66 

06 

17.0 

K5 

19.60 

1.29 

2.31 

0.76 

0.34 

14.59 





14.00 

14.08 



8.3 




IHq 764 

21 

42 

29.57 

+66 

05 

20.9 

G8-K2 

19.72 

0.96 

1.80 

0.50 

0.82 

15.25 









7.5 




IHq 765 

21 

42 

32.40 

+66 

04 

59.1 

G5-K0 

18.57 

0.86 

1.63 

0.56 

0.08 

15.01 





14.52 

15.12 



11.5 

ecr 



IHq 766 

21 

42 

34.72 

+66 

05 

18.6 

M3V 

17.80 

1.17 

2.54 

0.64 

0.15 

13.13 

12.87 

12.80 

12.82 

13.06 

13.22 

13.15 



4.3 

ecr 

S3-X4 


MMN 2 

21 

42 

38.80 

+66 

06 

35.8 


20.68 

2.22 

3.60 

1.32 

0.97 

12.51 

11.01 

10.59 

10.42 

10.19 

10.96 

10.49 



437 


S3-U939, HH 242 

wk cont 

S3-X41 

21 

42 

38.91 

+66 

07 

08.7 


18.37 

1.67 

3.20 

1.09 

0.44 

11.08 





10.42 

10.33 







IHq 767 

21 

42 

40.23 

+66 

13 

28.7 





0.94 

0.58 

12.51 





11.74 

11.05 

8.92 

6.70 

206 



wk cont 

S3-X52 

21 

42 

40.34 

+66 

10 

07.2 

A1 

12.44 

0.33 

0.73 

0.30 

0.13 

10.58 

10.52 

10.52 

10.27 

9.73 

10.43 

10.37 





SVS 2 


IHq 768 

21 

42 

40.49 

+66 

09 

51.6 


19.95 

0.85 

1.58 



14.87 

13.12 

12.33 

11.54 

10.48 

13.13 

11.92 



86.0 


S3-U1612 


MMN 3 

21 

42 

41.92 

+66 

09 

24.5 


21.19 

1.82 

3.78 

0.84 

0.35 

13.99 

13.56 

13.29 

13.08 

12.30 

13.66 

13.23 



22: 


S3-U1522 

wk cont 

IHq 769 

21 

42 

44.22 

+66 

10 

07.0 

M4.5V 

21.28 

1.82 

3.82 

0.83 

0.24 

14.30 





13.80 

13.29 



8.4 

0.8 



S3-X40 

21 

42 

45.32 

+66 

07 

04.4 


20.91 

1.63 

3.42 

1.02 

0.59 

13.39 

13.13 

13.05 

12.35 


11.39 

10.90 







IHq 770 

21 

42 

45.51 

+66 

06 

30.6 

GO; 

19.65 

1.18 

2.26 

0.83 

0.31 

14.41 









5.0 

abs 



IHq 771 

21 

42 

46.08 

+66 

05 

56.2 

M2V 

18.17 

1.31 

2.60 

0.77 

0.27 

12.75 

12.58 

12.60 

12.57 

13.15 

12.56 

12.38 



5.3 

abs 

S3-X31 


BD+65°1636 

21 

42 

46.09 

+66 

05 

13.8 

B8V 

10.86 

0.21 

0.39 

0.16 

0.07 

9.83 

9.85 

9.84 

9.87 

9.91 

9.66 

9.62 





S3-X3 


IHq 772 

21 

42 

46.12 

+66 

13 

45.8 





0.72 

0.41 

13.55 





13.29 

12.92 

10.88 


50.1 




IHq 773 

21 

42 

46.17 

+66 

06 

56.6 

MOV 

19.65 

1.45 

2.86 

0.77 

0.27 

12.75 





12.56 

12.38 



9.5 

0.2 



IHq 774 

21 

42 

46.55 

+66 

06 

33.4 

K3 

19.42 

1.24 

2.16 

0.70 

0.56 

14.35 









4.2 

abs 



IHq 775 

21 

42 

46.59 

+66 

06 

22.3 

M5V 

21.63 

1.30 

3.33 

0.74 

0.40 

14.98 









8: 

flat 



IHq 776 

21 

42 

46.87 

+66 

06 

57.4 

MO 

19.65 

1.45 

2.86 

1.02 

0.51 

12.09 

11.24 

10.91 

10.55 

9.74 





25.9 

abs 

S3-X39 


IHq 777 

21 

42 

47.05 

+66 

04 

57.8 


15.54 

1.02 

2.00 

0.77 

0.28 

11.05 

10.79 

10.74 

10.54 

10.00 

10.67 

10.37 



2.5 


S3-X2 


IHq 778 

21 

42 

47.46 

+66 

07 

03.4 

M2V 

21.66 

1.72 

3.84 

1.02 

0.52 

13.47 









20.8 

0.7 



IHq 779 

21 

42 

47.90 

+66 

06 

53.0 


17.91 

1.28 

2.53 

1.04 

0.33 

12.09 

11.69 

11.71 

10.89 






4.4 


S3-X23 


IHq 780 

21 

42 

48.13 

+66 

07 

43.2 

K5: 

21.42 

1.76 

3.22 

1.01 

0.37 

14.31 









4.9 

abs 



IHq 781 

21 

42 

49.87 

+66 

05 

42.6 


20.02 

1.36 

3.29 

0.64 

0.29 

13.72 

13.53 

13.39 







4: 


S3-X8 

em? 

IHq 782 

21 

42 

49.93 

+66 

05 

54.6 

M4 

21.61 

2.10 

4.08 

0.74 

0.45 

14.09 









18.3 

abs 



BD+65°1637 

21 

42 

50.18 

+66 

06 

35.1 

B3e 

10.53 

0.24 

0.51 

0.24 

0.26 

8.47 

8.13 

7.90 

7.57 

7.09 

7.90 

7.51 

6.17 

0.34 

19.6 


S3-X13 


IHq 783 

21 

42 

50.92 

+66 

06 

03.6 

M2V 

17.73 

1.21 

2.38 

0.85 

0.24 

12.55 

12.42 

12.37 

12.24 






5.3 

abs 

S3-X11 


IHq 784 

21 

42 

50.99 

+66 

03 

59.2 


19.96 

1.01 

1.88 








15.22 

14.07 



8: 



wk cont 

MMN 5 

21 

42 

51.42 

+66 

05 

56.2 


21.82 

1.56 

3.03 

1.07 

0.56 

13.57 

12.62 

12.09 

11.52 

10.53 





280 


S3-U815 

wk cont 

S3-X35 

21 

42 

51.96 

+66 

06 

33.4 


18.06 

0.99 

2.39 

0.73 

0.18 

13.13 

12.63 

12.68 











IHq 785 

21 

42 

52.30 

+66 

05 

35.3 

MIV 

18.98 

1.30 

2.70 

0.86 

0.32 

13.19 









3.7 

abs 



MMN 6 

21 

42 

52.61 

+66 

06 

57.2 

MIV 

18.95 

1.48 

2.94 

1.17 

0.85 

11.80 

10.81 

10.24 

9.77 

9.15 





69.2 

2.5 

S3-X25 


MMN 7 

21 

42 

53.14 

+66 

07 

14.8 


20.58 

1.12 

1.94 

1.42 

1.23 

12.93 

11.61 

10.93 

10.41 

9.57 





201 

24.6 

S3-U1085 

wk cont 

IHq 786 

21 

42 

53.21 

+66 

07 

20.8 


19.30 

1.42 

2.71 

1.15 

0.51 

12.68 

11.62 

11.29 

10.82 

9.92 





4.2 


S3-U1109 


MMN 8 

21 

42 

53.45 

+66 

09 

19.6 


21.87 

1.71 

2.77 

1.49 

0.92 

14.58 

12.91 

12.36 

11.87 

11.23 

13.25 

12.47 



234 


S3-U1504 

wk cont 

MMN 9 

21 

42 

53.49 

+66 

08 

05.3 

M2V 

17.55 

1.26 

2.50 

0.85 

0.24 

12.12 

11.85 

11.78 

11.22 


11.22 

11.16 



2.6 

abs 

S3-X29 


S3-X45 

21 

42 

54.08 

+66 

08 

14.9 

G8 

19.80 

1.57 

2.98 

1.17 

0.43 

12.82 

12.52 

12.55 

12.55 

12.29 









IHq 787 

21 

42 

54.71 

+66 

06 

35.6 

MOV 

19.00 

1.43 

2.59 

1.00 

0.62 

12.56 

11.66 

11.37 

11.02 

10.05 





5.4 

abs 

S3-X36 


MMN 10 

21 

42 

54.80 

+66 

06 

12.7 


19.01 

1.38 

2.59 

0.94 

0.35 

12.91 

12.58 








em? 


S3-X32 

wk cont 

IHq 788 

21 

42 

54.87 

+66 

06 

31.4 

MIV 

19.42 

1.53 

2.90 

1.05 

0.27 

13.42 









4.9 

abs 



IHq 789 

21 

42 

54.89 

+66 

07 

21.3 

M3V 

20.00 

1.71 

3.40 

1.05 

0.35 

12.84 

12.50 

12.44 

11.96 






5.8 

abs 

S3-X26 


S3-X37 

21 

42 

55.72 

+66 

06 

45.0 


21.16 

1.61 

3.64 

1.00 

0.31 

14.27 

13.62 

13.37 











IHq 790 

21 

42 

55.76 

+66 

05 

42.9 

M3.5V 

20.41 

0.74 

2.81 

0.85 

0.40 

13.97 

13.43 

13.38 







5.7 

ecr 

S3-X48 


IHq 791 

21 

42 

55.86 

+66 

07 

21.1 

MIV 

20.70 

1.69 

3.28 

1.31 

0.82 

12.47 

11.16 

10.54 

9.98 

9.02 





28.1 

8.6 

S3-U1107 


MMN 11 

21 

42 

56.25 

+66 

06 

02.1 


15.75 

0.95 

1.85 

0.77 

0.25 

11.41 

11.28 

11.13 

10.78 






em? 


S3-X10 


S3-X18 

21 

42 

56.77 

+66 

06 

37.2 


17.93 

1.31 

2.57 

0.98 

0.48 

12.21 

11.47 

11.21 

10.44 










IHq 792 

21 

42 

57.15 

+66 

06 

34.9 

K5 

17.18 

1.11 

2.11 

0.92 

0.22 

12.53 


12.00 







1.6 

abs 

S3-X15 


IHq 793 

21 

42 

57.40 

+66 

07 

15.1 

M3V 

19.72 

1.69 

3.38 

0.81 

0.36 

12.98 









3.2 

abs 



MMN 12 

21 

42 

58.10 

+66 

07 

39.3 

M5V 

20.49 

1.94 

3.93 

0.77 

0.50 

13.03 









13.6 

flat 



IHq 794 

21 

42 

58.18 

+66 

05 

40.5 

M2V 

18.64 

1.09 

2.27 

0.80 

0.13 

13.20 









8.9 

4.8 

S3-X47 


S3-X43 

21 

42 

58.34 

+66 

07 

26.3 


22.68 

2.53 

4.93 

1.69 

0.71 

11.91 

11.44 

11.39 

10.93 










IHq 795 

21 

42 

58.36 

+66 

05 

27.3 

M2.5V 

18.88 

1.30 

2.75 

0.93 

0.39 

12.64 

11.79 

11.39 







18.6 

0.5 

S3-X6 


BD+65°1638 

21 

42 

58.60 

+66 

06 

10.3 

B3V 

10.42 

0.33 

0.69 

0.32 

0.15 

8.51 

8.50 

8.48 

8.42 

8.00 







S3-X13, SVS 8 


S3-X17 

21 

42 

58.77 

+66 

06 

36.8 


17.84 

1.19 

2.46 

0.93 

0.22 

12.34 

11.98 

11.85 











IHq 796 

21 

42 

59.41 

+66 

11 

12.3 





0.76 

0.53 

14.77 









50: 



wk cont 

HBC 731 

21 

42 

59.61 

+66 

04 

33.8 

M1.5V 

18.69 

1.58 

3.08 

1.26 

0.68 

10.95 

10.18 

9.73 

9.39 

8.91 

10.20 

9.59 

7.58 

5.60 

54.9 

4.5 

S3-X1 


IHq 798 

21 

42 

59.97 

+66 

01 

01.0 





0.84 

0.20 

13.67 





13.60 

13.47 



em? 




IHq 799 

21 

42 

59.98 

+66 

06 

42.4 

M3V 

19.69 

1.47 

3.26 

0.86 

0.34 

13.15 

12.73 

13.05 







7.0 

ecr 

S3-X19 


HBC 732 

21 

43 

00.00 

+66 

11 

27.9 

M2 

16.35 

1.01 

2.11 

1.02 

0.68 

11.01 





10.06 

9.13 

6.58 

4.10 

27.7 


MMN 13, V350 Cep 


MMN 14 

21 

43 

00.23 

+66 

06 

47.4 

MOV 

20.02 

1.53 

2.73 

1.38 

0.69 

12.09 





9.53 

9.33 



92.7 

9.3 



SVS 13 

21 

43 

01.71 

+66 

07 

08.9 

B5 

14.57 

0.86 

1.78 

0.68 

0.31 

10.25 

9.19 

8.81 



7.64 

7.14 



em 

ecr 

S3-X51 


IHq 800 

21 

43 

01.88 

+66 

06 

44.7 

K3-5V 

17.24 

1.18 

2.22 

0.89 

0.40 

12.02 

11.12 

10.80 

11.15 






3.8 

abs 

S3-X20 




























Table 1 —Continued 


Identifier®' 

(J2000) 

5 

(J2000) 

SpT 

Vb 

V - Rpb 

V - 1+ 

J - 

H - Ks<= 

K'S' 

[3.6]<i 

[4.5]<1 

[5.81‘l 

[S.O]"* 

wl^ 


u;3® 

tt)4® 

(A) 

W(8542* 

(A) 

Other Identifiers 

Notes^ 

MMN 15 

21 

43 

02.46 

+ 66 

07 

03.9 


19.64 

1.34 

2.55 

0.90 

0.85 

12.26 









9: 




S3-X50 

21 

43 

03.01 

+ 66 

06 

55.9 








12.79 

12.59 











IHq 801 

21 

43 

03.20 

+ 66 

11 

15.0 





1.16 

1.12 

14.34 





11.70 

9.90 

6.01 

3.33 

111 


GGD 33A 


IHq 802 

21 

43 

03.43 

+ 66 

05 

26.4 

M2 

17.86 

1.24 

2.49 

0.98 

0.53 

12.01 









59.3 

4.5 



IHa 803 

21 

43 

04.71 

+ 66 

00 

30.5 





1.16 

0.62 

13.51 





12.93 

12.31 

10.68 


28: 



wk cont 

IHq 804 

21 

43 

04.95 

+ 66 

06 

53.6 

K2V 

17.90 

1.42 

2.58 

1.08 

0.72 

11.26 









6.3 

abs 

S3-X22 


IHq; 805 

21 

43 

05.09 

+ 66 

09 

29.5 

M7 

22.45 

1.93 

4.05 

0.64 

0.26 

14.90 

14.23 

14.07 

13.70 

13.11 

14.62 

14.33 



14.9 

abs 

S3-U1542 


LkHo 234 

21 

43 

06.82 

+ 66 

06 

54.2 

B8 

12.48 

0.66 

1.34 

1.33 

1.12 

7.08 

5.79 

5.15 

4.45 

3.33 

4.69 

2.90 

1.32 

-1.75 

56.6 

9.3 

MMN 16 


IHq 806 

21 

43 

11.41 

+ 66 

12 

55.5 





1.02 

0.60 

12.70 





11.89 

11.26 

9.43 

7.32 

166 




MMN 17 

21 

43 

11.61 

+ 66 

09 

11.4 

MIV 

16.30 

1.08 

2.16 

0.80 

0.31 

11.49 

11.04 

10.73 

10.38 

9.70 

11.16 

10.79 

8.82 

7.03 

4.8 

abs 

S3-X30 


IHq 807 

21 

43 

12.26 

+ 66 

06 

05.8 

M5V 

21.65 

1.48 

3.80 

0.74 

0.40 

13.95 

13.36 

13.05 

12.48 






7.0 

flat 

S3-U849 


IHq 808 

21 

43 

12.35 

+ 66 

09 

55.4 

M4.5 

19.71 

1.52 

3.08 

0.91 

0.31 

13.59 

13.00 

12.74 

12.39 

11.85 

13.15 

12.73 

11.18 


17.4 

ecr 

S2-U1640 


IHq 809 

21 

43 

14.36 

+ 66 

08 

59.2 

M4V 

22.48 

1.62 

4.21 

0.73 

0.37 

14.93 

14.23 

13.89 

13.59 

12.98 





36.9 

abs 

S3-U1433 


IHq 810 

21 

43 

14.43 

+ 66 

07 

34.8 

M3 

20.50 

1.41 

2.95 

1.18 

0.49 

13.84 





12.89 

10.98 



2.9 

flat 



IHq 811 

21 

43 

15.28 

+ 66 

07 

57.1 

M2.5V 

21.31 

1.66 

3.38 

1.67 

0.75 

12.59 

11.57 

11.03 

10.41 

9.41 





6.2 

flat 

S3-X28 


MMN 19 

21 

43 

16.83 

+ 66 

05 

48.6 

M2V 

18.31 

1.29 

2.48 

0.88 

0.18 

13.12 

12.74 

12.57 

12.45 

11.99 





15.8 

abs 

S3-X9 


IHq 812 

21 

43 

18.06 

+ 66 

05 

35.2 

M3V 

19.74 

1.47 

3.06 

0.77 

0.29 

13.70 

13.39 

13.27 

13.29 

13.12 





3.6 

flat 

S3-X7 


IHq 813 

21 

43 

19.37 

+ 66 

07 

21.5 

M4.5V 

20.70 

1.67 

3.50 

0.77 

0.40 

13.89 

13.90 

13.72 

13.62 

13.43 





17.5 

0.8 

S3-X42 


IHq 814 

21 

43 

21.07 

+ 66 

06 

22.8 


17.83 

1.23 

2.56 

0.66 

0.15 

12.75 

12.44 


12.32 

12.42 





2; 

abs 

S3-X14 

binary 

S2-X2 

21 

43 

24.89 

+ 66 

07 

34.1 


23.25 

1.44 

3.00 




17.12 

16.15 











IHq 815 

21 

43 

26.95 

+ 66 

09 

36.5 

MOV 

16.49 

0.87 

1.67 

0.71 

0.18 

12.72 

12.67 

12.66 

12.60 

12.55 

12.73 

12.74 



1.6 

ecr 

S2-X4 


IHq 816 

21 

43 

28.10 

+ 66 

00 

57.2 





0.95 

0.33 

12.55 





12.25 

11.91 

10.77 


15.0 




S3-X56 

21 

43 

29.28 

+ 66 

03 

31.5 

G6 

11.17 

0.31 

0.58 

0.29 

0.04 

9.69 

9.69 

9.69 

9.64 

9.64 

9.64 

9.65 

10.13 

9.41 



SVS 16 


MMN 20 

21 

43 

31.82 

+ 66 

08 

50.6 


18.21 

1.20 

2.34 

1.04 

0.60 

12.30 

11.59 

11.29 

11.02 

10.31 

11.88 

11.47 



21.1 


S2-X1 


IHq 817 

21 

43 

36.87 

+ 66 

07 

52.6 

K5: 

19.33 

1.16 

2.19 

0.70 

0.45 

14.40 





14.23 

14.19 



4.4 

abs 



MMN 22 

21 

43 

43.44 

+ 66 

07 

30.8 


18.57 

1.38 

2.70 

1.05 

0.46 

12.15 

11.43 

11.05 

10.63 

10.06 

11.65 

11.16 

9.60 


em? 


S2-X5 


IHq 818 

21 

43 

43.71 

+ 66 

08 

22.3 

M3 

19.61 

1.45 

2.95 

0.90 

0.25 

13.66 









7.3 

flat 



IHq 819 

21 

44 

05.37 

+ 66 

05 

53.1 

K5 

18.70 

1.09 

2.15 

1.54 

0.93 

11.22 

9.83 

9.33 

8.92 

8.16 

10.09 

9.38 

7.32 

4.94 

60.6 

5.3 

S2-U820 


IHq 820 

21 

44 

06.34 

+ 66 

04 

23.1 


18.71 

1.36 

2.62 

1.13 

0.69 

11.38 





10.51 

10.00 

7.96 

5.68 

2.7 




IHq 821 

21 

44 

07.82 

+ 66 

04 

33.2 

M3.5 

19.57 

1.59 

3.16 

0.94 

0.31 

13.21 





12.70 

12.26 



25.4 

ecr 




^IHq number unless previously identified by Magakian et al. (2004), Herbig (1957), or the Herbig Bell Catalog (HBC). 

^Optical photometry from the UH 2.2 m or the KPNO 0.9 m. 

‘"Near-infrared photometry from the 2MASS Point Source Catalog. 

Spitzer IRAC photometry from Stelzer & Scholz (2009). 

®WISE photometry from the AllWISE Source Catalog. 

^Positive values indicate emission. 

^Positive values indicate emission; abs - Ca II A8542 in absorption; ecr - emission core reversal is evident; flat - flat continuum noted at A8542 
^wk cont: weak continuum on WFGS images. 













Table 2 

Photometry for Infrared Sources without X-ray or Ha Emission^ 


Source 

(J2000) 

S 

(J2000) 

Rc - 


J - 

H - 


[3.6]<^ 

[4.5]'^ 

[5.8]'* 

[8.0]<i 


-u;2® 

■u>3® 

u;4® 

Other Identifiers 

S3-U1178 

21 41 55.30 

+ 66 

07 41.5 






15.86 

14.90 

13.98 

12.78 

16.10 

15.12 

12.15 

8.99 


S3-U1367 

21 42 17.67 

+ 66 

08 40.3 



1.15 

0.39 

14.63 

13.68 

13.23 

12.91 

12.40 

14.06 

13.61 




S3-U1194 

21 42 42.42 

+ 66 

07 45.2 

2.42 

18.62 

1.35 

0.78 

13.45 

12.56 

12.21 

12.23 

11.68 

12.07 

11.80 




S3-U1521 

21 42 42.86 

+ 66 

09 24.0 






16.23 

15.25 

13.84 

12.79 






S3-U1246 

21 42 48.23 

+ 66 

08 00.6 

2.30 

18.16 

1.40 

0.81 

13.08 

11.69 

11.14 

10.62 

9.67 

11.23 

10.71 




S3-U722 

21 42 54.63 

+ 66 

05 20.3 

2.66 

18.50 

0.95 

0.69 

13.87 

12.51 

12.64 








S3-U546 

21 42 57.75 

+ 66 

04 23.5 

1.08 

18.00 


0.82 

13.66 

11.80 

10.94 

9.95 

8.54 

11.85 

10.74 

7.10 

5.07 

V392 Cep, RNO 138 

S3-U1780 

21 42 59.47 

+ 66 

10 35.8 




1.07 

14.61 

13.54 

13.00 

12.52 

12.05 

13.88 

13.09 




S3-U270 

21 42 59.82 

+ 66 

01 54.9 






16.83 

16.59 

15.91 

14.33 






S3-U419 

21 43 01.78 

+ 66 

03 24.4 






13.02 

10.41 

10.17 

9.08 

13.17 

9.86 

6.03 

0.36 

FIRS 2 

S3-U500 

21 43 02.01 

+ 66 

04 02.7 






13.64 

12.98 

12.18 

11.29 

13.37 

12.32 




S3-U1611 

21 43 02.62 

+ 66 

09 50.7 

2.60 

18.08 

0.92 

0.59 

14.20 

13.54 

13.16 

12.65 

11.93 

13.92 

13.32 




S3-U1294 

21 43 02.89 

+ 66 

08 14.2 



0.92 

0.41 

15.01 

14.15 

13.83 

13.58 

13.43 






S3-U822 

21 43 04.38 

+ 66 

05 56.4 


18.55 

0.58 

0.65 

14.47 


12.55 








S3-U1242 

21 43 05.92 

+ 66 

07 58.5 






13.78 

12.70 

12.10 

11.65 






S3-U968 

21 43 06.96 

+ 66 

06 41.7 




2.34 

10.89 

7.54 

6.59 

5.67 

5.06 

8.38 

6.75 

4.15 

0.70 

J21430696-I-660641.7 

S3-U1103 

21 43 07.84 

+ 66 

07 18.5 

2.42 

19.38 

1.87 

0.91 

13.94 

13.04 

12.57 

11.76 

10.88 






S3-U550 

21 43 11.17 

+ 66 

04 25.6 






13.91 

13.42 

12.98 

12.33 

13.83 

13.11 




S3-U1026 

21 43 12.42 

+ 66 

06 55.9 






14.16 

13.73 

13.25 

12.23 






S3-U1211 

21 43 14.17 

+ 66 

07 46.5 





14.34 

10.55 

9.34 

8.41 

7.51 






S3-U1169 

21 43 14.83 

+ 66 

07 37.5 





14.77 

12.93 

10.72 

10.44 

10.22 






S3-U1350 

21 43 24.13 

+ 66 

08 31.5 




1.43 

14.19 

12.86 

11.28 

10.05 

8.96 

13.07 

11.04 

8.85 

4.64 

GGD 34C 

S3-U1059 

21 43 24.90 

+ 66 

07 04.7 

1.58 

19.72 




14.23 

13.33 

12.59 

11.56 






S2-U1313 

21 43 26.64 

+ 66 

08 20.5 






14.66 

14.19 

13.79 

13.24 

14.97 

13.92 




S3-U821 

21 43 29.32 

+ 66 

05 55.7 

2.70 

18.46 

0.92 

0.57 

14.39 

13.47 

13.04 

12.52 

11.02 

13.71 

13.07 




S2-U1083 

21 43 29.92 

+ 66 

07 09.2 


21.79 




14.67 

14.24 

13.99 

13.38 






S2-U804 

21 43 33.12 

+ 66 

05 49.0 






16.69 

15.83 

14.92 

13.37 






S2-U2219 

21 43 38.59 

+ 66 

12 30.6 






13.58 

13.34 

12.80 

11.99 

13.42 

13.06 




S2-U1660 

21 43 47.04 

+ 66 

10 01.2 

1.11 

17.04 

0.70 

0.51 

14.19 

14.05 

14.03 

13.72 

13.00 

13.84 

13.82 




S2-U613 

21 43 48.70 

+ 66 

04 46.2 



1.23 

0.57 

13.08 

12.40 

12.15 

11.83 

10.87 

12.33 

11.88 

9.65 



S2-U1713 

21 43 49.35 

-^66 

10 12.9 

1.24 

17.79 

0.79 

0.21 

14.86 

14.53 

14.47 

14.29 

13.36 

14.33 

14.53 





^Young Stellar Objects identified by Stelzer Scholz (2009) based upon infrared excess emission. 
^Optical photometry from Keck LRIS imaging, UH 2.2 m or the KPNO 0.9 m. 

^Near-infrared photometry from the 2MASS Point Source Catalog. 

^Spitzer mid-infrared photometry from Stelzer &; Scholz (2009). 

®Mid-infrared photometry from the AllWISE Source Catalog. 













Table 3 

Photometry for Spectroscopically Classified Stars without X-ray or Ha Emission 


Source^ 

(J2000) 

S 

(J2000) 

SpT 


V - 

Rc - ^c*’ 

J - 

H - Kg' 



tu2^ 

•u;3‘^ 

1114*^ 

Comments 

1 

21 42 22.75 

+ 66 07 52.9 

G3-5 

16.65 

0.82 

0.75 

0.54 

0.22 

13.25 

13.06 

13.14 




2 

21 42 22.85 

+ 66 06 56.3 

G-Kl 

20.57 

1.34 

1.11 

0.80 

0.55 

14.90 


14.06 



IR excess 

3 

21 42 23.07 

+ 66 06 42.5 

Kl-3 

18.86 

1.04 

0.87 

0.66 

0.13 

14.59 






4 

21 42 23.31 

+ 66 08 47.4 

MO 

17.45 

1.01 

0.93 

0.68 

0.21 

13.45 

13.35 

13.35 




5 

21 42 24.74 

+ 66 06 21.4 

P9-G5 

17.59 

1.23 

1.14 

0.88 

0.32 

12.22 

11.72 

11.88 




6 

21 42 25.19 

+ 66 10 09.0 

F7-G5 

18.14 

0.96 

0.91 

0.42 

0.32 

14.15 

13.99 

13.92 




7 

21 42 26.84 

+ 66 07 42.5 

G9-K0 

10.59 

0.57 

0.50 

0.55 

0.07 

7.96 

7.84 

7.93 

7.68 


SVS 4 

8 

21 42 28.35 

+ 66 04 08.7 

P5-9 

18.41 

0.99 

0.95 

0.62 

0.22 

14.18 

13.43 

13.66 



IR excess 

9 

21 42 31.22 

+ 66 08 04.8 

G9-K1 

19.17 

1.01 

0.91 

0.52 

0.29 

14.98 

14.55 





10 

21 42 31.86 

+ 66 07 08.7 

G6-8 

17.29 

0.99 

0.95 

0.61 

0.18 

13.11 

12.53 

12.71 




11 

21 42 32.57 

+ 66 10 24.9 

G8-K0 

18.54 

1.02 

0.96 

0.69 

0.16 

14.26 

14.07 

14.23 




12 

21 42 33.75 

+ 66 08 02.3 

M2.5 

19.93 

1.69 

1.55 

0.93 

0.40 

13.37 

12.98 

13.03 




13 

21 42 34.28 

+ 66 03 12.6 

G-K4 

18.37 

1.20 

0.92 

0.79 

0.23 

13.73 

13.42 

13.59 




14 

21 42 35.03 

+ 66 04 09.6 

G5-K0 

18.36 

0.82 

0.77 

0.52 

0.23 

14.90 

14.10 

14.99 




15 

21 42 35.69 

+ 66 05 49.3 

K3 

17.44 

1.05 

0.95 

0.74 

0.32 

12.80 

12.77 

13.02 




16 

21 42 38.58 

+ 66 09 41.4 

M2.5 

20.27 

1.35 

1.32 

0.50 

0.70 

15.10 

14.39 

14.26 



IR excess 

17 

21 42 39.05 

+ 66 07 09.9 

K5 

18.37 

1.67 

1.53 

1.09 

0.44 

11.08 

10.46 

10.39 



IR excess 

18 

21 42 41.26 

+ 66 03 33.7 

G2-5 

16.02 

0.77 

0.74 

0.52 

0.12 

12.80 

12.66 

12.73 

12.27 



19 

21 42 44.96 

+ 66 07 48.0 

G2-K 

22.03 

2.06 

1.57 

1.03 

0.48 

14.06 






20 

21 42 45.98 

+ 66 08 15.6 

K5 

20.49 

1.68 

1.49 

1.03 

0.52 

13.44 

12.57 

12.49 



IR excess 

21 

21 42 46.02 

+ 66 10 20.3 

K5: 

22.69 

2.18 

2.21 

1.11 

0.65 

13.34 

12.79 

12.39 



IR excess 

22 

21 42 47.73 

+ 66 05 35.2 

G2-5 

17.86 

0.97 

0.86 

0.64 

0.08 

13.86 

12.84 

12.54 



IR excess 

23 

21 42 50.66 

+ 66 03 31.3 

G9-K2 

16.95 

0.73 

0.58 

0.61 

0.17 

13.99 

13.88 





24 

21 42 50.47 

+ 66 08 32.7 

K2-5 

21.26 

1.53 

1.37 

1.32 

0.31 

14.74 






25 

21 42 51.16 

+ 66 05 45.3 

P2 

17.29 

0.59 

0.58 

0.26 

-0.03 

14.90 






26 

21 42 51.94 

+ 66 09 44.9 

Ml 

18.05 

1.03 

0.91 

0.70 

0.23 

13.93 

13.97 





27 

21 42 53.97 

+ 66 08 15.1 

G5-K0 

19.80 

1.57 

1.41 

1.17 

0.43 

12.82 






28 

21 42 55.18 

+ 66 06 25.4 

K5 

15.83 

0.84 

0.68 

0.68 

0.15 

12.53 






29 

21 42 58.04 

+ 66 05 54.7 

G3-5 

14.21 

0.41 

0.40 

0.31 

0.20 

12.26 






30 

21 43 00.20 

+ 66 05 46.1 

G8 

17.47 

0.76 

0.65 

0.84 

0.01 

13.87 






31 

21 43 05.32 

+ 66 04 52.0 

Ml 

19.57 

1.09 

1.11 

0.74 

0.30 

14.93 






32 

21 43 05.05 

+ 66 10 12.4 

Ml 

20.97 

1.37 

1.51 

1.02 

0.48 

14.14 

14.07 

14.03 




33 

21 43 05.53 

+ 66 03 28.7 

K5 

16.93 

0.86 

0.66 

0.72 

0.07 

13.63 

13.41 

13.45 




34 

21 43 07.51 

+ 66 09 01.5 

M3 

20.84 

1.85 

1.72 

1.12 

0.36 

13.20 






35 

21 43 09.85 

+ 66 05 48.8 

M3 

22.52 

2.46 

1.93 

1.02 

0.47 

13.86 






36 

21 43 13.04 

+ 66 05 50.8 

M4.5 

20.49 

1.62 

1.84 

0.80 

0.39 

13.64 






37 

21 43 14.63 

+ 66 10 13.6 

K1 

18.46 

1.20 

1.13 

0.78 

0.29 

13.31 

13.16 

13.10 




38 

21 43 14.93 

+ 66 09 07.1 

MO 

17.99 

1.05 

0.92 

0.66 

0.26 

13.86 

13.73 

13.80 




39 

21 43 17.55 

+ 66 04 18.5 

GO-2 

20.29 

1.51 

1.42 

0.89 

0.47 

13.58 






40 

21 43 20.84 

+ 66 03 37.1 

Ml 

17.55 

1.04 

1.01 

0.68 

0.19 

13.45 

13.22 

13.19 




41 

21 43 21.66 

+ 66 02 46.1 

G3 

14.25 

0.44 

0.42 

0.40 

0.03 

12.40 

12.32 

12.34 




42 

21 43 31.12 

+ 66 07 43.7 

M2 

19.38 

1.18 

1.18 

0.61 

-0.04 

14.90 


14.44 




43 

21 43 31.18 

+ 66 07 24.4 

G8-K1 

14.53 

0.62 

0.52 

0.52 

0.12 

11.96 

11.86 

11.93 




44 

21 43 31.18 

+ 66 09 54.7 

G6-8 

15.70 

1.33 

1.24 

0.93 

0.32 

9.96 

9.74 

9.72 

9.85 


SVS 14 

45 

21 43 32.68 

+ 66 10 11.8 

G-Kl 

20.16 

1.49 

1.38 

1.19 

0.25 

13.55 

13.15 

13.15 




46 

21 43 40.01 

+ 66 03 31.9 

G4 

12.92 

0.38 

0.34 

0.33 

0.04 

11.13 

11.06 

11.09 

11.20 



47 

21 43 40.64 

+ 66 05 16.5 

M5 

21.50 

1.38 

1.80 

0.56 

0.76 

15.08 

15.14 

15.05 




48 

21 43 44.89 

+ 66 07 0.12 

A2 

15.29 

0.90 

0.90 

0.49 

0.21 

11.40 

11.17 

11.13 




49 

21 43 45.10 

+ 66 04 38.1 

G8 

15.57 

0.56 

0.46 

0.37 

0.15 

13.28 

13.18 

13.21 




50 

21 43 45.37 

+ 66 08 22.9 

Ml 

17.50 

1.09 

1.10 

0.58 

0.28 

13.22 

12.83 

12.62 



IR excess 

51 

21 43 47.03 

+ 66 10 01.7 

G5-K0 

19.37 

1.22 

1.11 

0.70 

0.51 

14.19 

13.84 

13.82 




52 

21 43 47.40 

+ 66 08 49.5 

M3 

20.01 

1.16 

1.26 

0.81 

0.48 

15.05 






53 

21 43 50.02 

+ 66 07 52.3 

M3.5 

16.84 

1.22 

1.46 

0.60 

0.23 

11.94 

11.68 

11.53 




54 

21 43 50.48 

+ 66 09 38.8 

K2 

18.72 

1.29 

1.22 

0.90 

0.30 

13.11 

12.69 

12.74 




55 

21 43 51.63 

+ 66 03 39.0 

KO-3 

19.95 

1.20 

1.12 

0.92 

0.19 

14.70 

14.66 

14.71 




56 

21 43 54.23 

+ 66 06 42.5 

K5: 

21.18 

1.66 

1.37 

1.14 

0.23 

14.37 

13.91 

13.91 




57 

21 43 54.30 

+ 66 09 10.4 

K5: 

20.25 

1.36 

1.20 

0.96 

0.45 

14.48 

14.12 

14.32 




58 

21 43 55.69 

+ 66 05 29.6 

K5: 

19.63 

1.46 

1.33 

0.98 

0.37 

13.29 

12.99 

12.97 




59 

21 43 55.58 

+ 66 08 16.9 

P7-G5 

18.54 

1.11 

1.03 

0.72 

0.19 

13.96 

13.66 

13.69 




60 

21 43 55.98 

+ 66 03 05.3 

G6 

15.68 

0.62 

0.53 

0.56 

0.10 

13.11 

13.04 

13.08 




61 

21 43 58.34 

+ 66 07 11.7 

K5 

21.26 

1.72 

1.29 

0.91 

0.45 

14.73 

14.55 

14.56 




62 

21 43 59.96 

+ 66 04 37.0 

K4 

16.90 

0.83 

0.65 

0.65 

-0.03 

13.83 

13.59 

13.67 




63 

21 44 00.52 

+ 66 05 17.3 

MO 

16.82 

1.00 

0.95 

0.74 

0.14 

12.85 

12.64 

12.63 




64 

21 44 00.95 

+ 66 07 24.7 

K1 

18.20 

1.40 

1.26 

0.97 

0.33 

12.31 

12.08 

12.11 















Table 3— 


Source®' 

(J2000) 

S 

(J2000) 

SpT 

v-b 

V - He*’ 

He - J-e*” 

65 

21 44 02.99 

4-66 04 56.2 

G-K 

19.20 

1.39 

1.23 

66 

21 44 04.32 

4-66 06 41.6 

K2 

16.74 

1.31 

1.20 

67 

21 44 05.12 

4-66 06 00.2 

P7 

17.04 

1.03 

0.96 

68 

21 44 05.41 

4-66 04 52.6 

M2 

18.65 

1.11 

1.09 

69 

21 44 06.59 

4-66 09 13.1 

G7-K0 

17.02 

1.14 

1.05 

70 

21 44 06.71 

4-66 08 10.0 

G5 

17.84 

0.96 

0.90 

71 

21 44 08.14 

4-66 06 04.2 

G7-K0 

18.51 

1.13 

1.02 

72 

21 44 08.64 

4-66 03 15.9 

G5-8 

19.47 

1.30 

1.18 


^Identifier 

^Optical photometry from the KPNO 0.9 m or the UH 2.2 m. 
‘'Near-infrared photometry from the 2MASS Point Source Catalog. 
‘^Mid-infrared photometry from the AllWISE Source Catalog. 


J - H' 


H - K 


Comments 


w2^ w3^ 


0.89 0.31 13.46 13.20 13.19 

0.90 0.27 11.23 11.02 11.04 

0.61 0.21 12.78 12.63 12.58 

0.55 0.16 14.50 14.20 14.15 

0.82 0.20 12.23 12.04 12.05 

0.56 0.19 14.01 13.87 13.85 

0.80 0.11 13.88 13.44 13.38 ... ... IR excess? 

0.86 0.42 13.91 13.70 13.70 









